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Abstract 

The factorization properties of the radiative decays R —> P 7 are analyzed at leading or¬ 
der in 1 /mb using the soft-collinear effective theory. It is shown that the decay amplitudes 
can be expressed in terms of a B ^ V form factor evaluated at = 0, light-cone distri¬ 
bution amplitudes of the B and V mesons, and calculable hard-scattering kernels. The 
renormalization-group equations in the effective theory are solved to resum perturbative 
logarithms of the different scales in the decay process. Phenomenological implications for 
the B —> K*j branching ratio, isospin asymmetry, and CP asymmetries are discussed, 
with particular emphasis on possible effects from physics beyond the Standard Model. 


1 Introduction 


In the Standard Model, the radiative b ^ s'y transition is suppressed and it is therefore 
a sensitive probe for the effects of New Physics. The total B X^'j decay rate can be 
calculated in an expansion about the heavy-quark limit using the operator product expansion 
(OPE). At leading order in the heavy-quark expansion, the total rate can be calculated in 
perturbation theory and it is therefore known rather precisely. 

However, the OPE is only valid for sufficiently inclusive observables. It cannot be used 
if the photon energy in the inclusive process is restricted to the endpoint region, much less 
to analyze the exclusive decay B —> K*^. Restricted to the region of large photon energy, 
the b ^ s transition involves non-perturbative strong-interaction physics, even in the heavy- 
quark limit. The factorization analysis retains predictive power by organizing these non- 
perturbative contributions in a universal and process-independent manner. An efficient way 
to study these decays in the heavy-quark expansion is to use soft-collinear effective theory 
(SCET) [1, 2, 3, 4, 5, 6]. In this approach, the relevant momentum regions in the pertinent 
Feynman diagrams are represented by helds in the effective theory, and the expansion of the 
diagrams in momentum space translates into a derivative expansion of the effective Lagrangian. 
The use of a Lagrangian makes the structure of the interactions and the resulting factorization 
properties of the amplitudes more transparent and allows identihcation of the remaining non- 
perturbative parts of a given process at the operator level. It provides a simple way to resum 
large perturbative logarithms associated with the different scales in the problem, by solving the 
renormalization-group (RG) equations obeyed by the effective-theory operators. For inclusive 
decay distributions in the endpoint region, SCET has been used to obtain a factorization 
theorem for the decay at next-to-leading order in the heavy-quark expansion [7, 8, 9], an 
analysis which seems prohibitively difficult on a purely diagrammatic level. 

In the present paper, we use SCET to analyze the exclusive decay B or more 

generally B ^ V'y, where R is a light vector meson. We demonstrate that, at leading order 
in 1/mi, and to all orders in as, the matrix elements of the operators Qi in the effective weak 
Hamiltonian governing the decay obey the generalized factorization formula 

—M^) du4,vJu)T‘\u,u). (1) 

Jo ^ Jo 

The quantities and are perturbatively calculable functions, appearing as Wilson co¬ 
efficients of effective-theory operators. jg form factor evaluated at maximum recoil 

= 0), and 0^, are light-cone distribution amplitudes (LCDAs) for the heavy and light 
mesons, respectively. 

As is manifest from the factorization theorem (1), the B ^ V form factors at large recoil 
energy are an important ingredient in the analysis of rare exclusive B ^ Vj decays. These 
form factors have been analyzed in the effective theory in [10, 11, 12]. It was found that 
the form factors in this energy regime contain a non-factorizable piece, which however is 
independent of the Dirac structure of the current in the heavy-quark limit. A single function 
Cv±{E) then suffices to describe all B ^ V± form factors up to factorizable corrections [13, 14]. 
The proof of the factorization theorem (1) is achieved after showing that the non-factorizable 
piece of the R —Ry decay amplitude is given by the same function. We will show that 
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diagrams in which the photon is emitted from one of the current quarks have the same structure 
as those encountered in the study of heavy-to-light form factors. Using the SCET formalism 
and the results from the form-factor analysis, it is then straightforward to establish (1) for these 
contributions. Photon emission from the 5-meson spectator quark has a more complicated 
structure. Diagrams of this type develop singularities for momentum conhgurations where 
some of the quarks and gluons are collinear to the photon (instead of being collinear to the 
meson V). Such conhgurations do not appear in the form-factor analysis, and to describe 
them it is necessary to include additional collinear helds in the effective theory. We introduce 
a counting scheme to systematically list all operators that may contribute at a given order in 
the power counting, and show that at leading power the operator matrix elements obey (1). 
The matching of QCD onto the effective theory is performed in two steps: at the hard scale 
fi ~ mb the operators are matched onto an intermediate theory, SCETi. The matching of 
SCETj onto the hnal effective theory, SCETn, is then performed at the hard-collinear scale 
H ~ Amt, where A is a typical hadronic scale. The two-step matching procedure makes it 
simpler to identify the operators appearing in SCETn and to separate the part proportional 
to the form factor from the remainder. Such a two-step matching is also required in order 
to resum large logarithms in perturbative expansions involving both the hard and the hard- 
collinear scales. The term in (1) involves only the hard scale, so that large logarithms 
are avoided by taking the scale jj, ~ mb- The term however involves both the hard and 
hard-collinear scales. In this case no scale choice is possible that avoids all large logarithms, 
and we perform the necessary resummation for this term. 

Operators describing spectator emission in radiative B decays appear at leading power in 
the effective theory, but in the Standard Model, the corresponding matrix elements between 
pseudoscalar B mesons and transversely-polarized vector mesons vanish. Such operators would 
contribute at leading power to the radiative decay B* P'y of B* into light pseudoscalar 
mesons. While this decay mode is not of phenomenological importance, it is interesting to 
note that, as we will show, the formula (1) extends without essential modihcation to the 
general case of radiative B or B* decay to a flavor non-singlet light meson M: 

r°° d( ; P 

/ —M^) du4)M{u)Tl\uj,u). ( 2 ) 

do ^ do 

Note that the form factors and heavy-quark LCDAs for the B and B* mesons are related in 
the heavy-quark limit. In the presence of New Physics operators with a chirality structure 
different from those of the Standard Model, the 5 —> Uy decay amplitude receives leading- 
power contributions associated with spectator photon emission [15], with both left- and right- 
circular photon polarization. We consider a general class of such operators and show that they 
also obey the factorization formula (1). Even at leading power in 1/mb-, these operators break 
isospin symmetry, and they can give rise to a non-vanishing time-dependent CP asymmetry. 
Measurements of these asymmetries in exclusive radiative B decays can thus provide useful 
constraints on the Wilson coefficients of the associated New Physics effective operators. For 
completeness, we consider also the case of flavor-singlet final-state mesons. Here we find new 
classes of operators not appearing in the form-factor analysis. These operators have vanishing 
matrix elements between pseudoscalar B mesons and transversely-polarized final-state vector 
mesons, but would in principle contribute to radiative B* decays to pseudoscalar final states. 


2 



One class of new operators is non-factorizable. Since these contributions cannot be related 
to form factors, we conclude that a factorization formula such as ( 2 ) does not hold for flavor 
singlet B* —P 7 decays. 

For the case of B ^ K*'y, we assess the impact of strange-quark mass effects on the 
predictions of the factorization theorem (1). We hrst demonstrate the theorem for the massless 
case, and then consider the perturbation caused by a small but hnite strange-quark mass, 

We argue that the leading corrections, linear in rris/A, simply contribute to the universal 
LCDA (pK*^ and non-factorizable form factor Cx*, breaking the SU(3) flavor symmetry that 
one obtains for the purely massless case. Possible corrections to the form of the factorization 
formula ( 1 ) itself could only appear starting at quadratic order, (mg/A)^. 

Factorization for the B ^ V'y decay process has received considerable attention in the 
literature. Extending the QCD factorization formalism [16] to the case of exclusive radiative 
decays, the formula (1) was proposed in [17, 18], where and were calculated through 
0{as)- Diagrams contributing one-loop matching corrections to were studied in [19], 
and (1) was shown to hold for such spectator interactions through one-loop order. Other 
studies include [20], and recent updates in [21, 22, 23]. Such explicit demonstrations give 
important insight into the structure of the decay process. However, an all orders proof of 
the validity of ( 1 ) is still lacking, an issue we address in this paper using the language of 
SOFT. In addition to this formal aspect of establishing factorization, the effective held-theory 
language has the advantage of systematically separating higher-order perturbative corrections 
that contribute to or SCET achieves this simplihcation by identifying the separate 
contributions on the operator level, a feature which also allows the resummation of large 
logarithms appearing in the perturbative kernels. A previous analysis of radiative B ^ Vj 
decay using SCET [24] identihed the SCETi operators corresponding to the contributions 
calculated in [17, 18]. However, to establish factorization one needs to construct a complete 
basis of effective-theory operators that can contribute to the process under consideration in 
the heavy-quark limit, an issue that was not addressed in [24]. Furthermore, this analysis 
suffers from an incomplete treatment of the low-energy theory, SCETn. As we will emphasize, 
a demonstration of factorization must deal with the soft-collinear messenger modes that can 
potentially spoil factorization [25]. In particular, the decoupling of soft gluons from hard- 
collinear helds in SCETi is not sufficient to ensure factorization; for example, the matrix 
element of the operator in [10, 24] is “factorizable” in this sense, but it cannot be written 
as a convergent convolution of a perturbative kernel with meson LCDAs [ 12 ]. 

The paper is organized as follows. In Section 2, we discuss the diagrammatic analysis of the 
decay. After identifying the necessary momentum regions we introduce the corresponding helds 
and set up the effective Lagrangian. In Section 3, we then hud the SCET operators needed to 
analyze the decay at leading power. This point needs special consideration, since the matching 
of SCETi onto SCETn involves inverse derivatives counting as negative powers of the expansion 
parameter. The most general operator at a given power is determined by dimensional analysis 
and longitudinal boost invariance. Aside from the contribution of messenger modes, which 
communicate between the soft and collinear sectors, this issue was addressed by Beneke and 
Feldmann in their analysis of heavy-to-light form factors [11]. We will extend their discussion 
to cover the more complicated case of photon emission from the spectator quark, involving 
two different types of collinear helds dehned with respect to opposite light-cone directions. 
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Using the same power-counting arguments we analyze the infrared messenger modes that 
can potentially spoil factorization in the matrix elements dehning the second term of (1). 
In Section 4, we match the effective weak Hamiltonian onto the list of operators derived 
in Section 3. The resulting B ^ V'j matrix elements can be written in the form of the 
factorization theorem (1). We show that the infrared messenger modes cannot contribute 
to the matrix elements dehning the second term of (1), thus demonstrating that the hard- 
scattering kernels are free of infrared divergences to all orders in perturbation theory. In 
Section 5 we consider the phenomenological implications of our analysis by computing the 
B — K*'y branching fraction, isospin asymmetry and CP asymmetry. We include the hrst 
complete treatment of the hard-scattering terms at leading order in RG improved perturbation 
theory. We discuss the phenomenological impact of the resummation of the leading single and 
double logarithms and of the inclusion of one-loop matching corrections to the hard-scattering 
kernel. In Section 6 we summarize our results and present our conclusions. 


2 Perturbative analysis oi B —> 


In the Standard Model, the effective weak Hamiltonian mediating havor-changing neutral 
current (FCNC) transitions of the type b ^ s has the form 


Hw — 


G 
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V2 ^ 

'' p=u^c 
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( 3 ) 
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(4) 


Here i and j are color indices. The effective weak Hamiltonian for b ^ d transitions is obtained 
by replacing s —d in the above expressions. Our sign conventions are such that the covariant 
derivative acting on a down-type quark is iD^ = id^ — \eA^ + 

Our task is to analyze the factorization properties of the matrix elements involving the 
above operators. The factorization theorem (1) holds trivially for the operator Qi, which 
directly maps onto the QCD tensor current. The goal of the present paper is to show that the 
matrix elements of the remaining operators can also be brought into this form. To analyze 
the factorization properties of these matrix elements, we use the reduction formula 



jd*y A) | T 
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(V(pv,y)7(Py,A)jQi(0) jB(pB)} + 


( 5 ) 


where the currents Jb and Jy have the quantum numbers of the B meson and the vector 
meson, respectively, with associated decay constants and fy\ The ellipsis stands for 
terms that do not have a pole at py = ruy with py > 0 and at p\ = m\ with p^ > 0. We 
then analyze the correlator on the left-hand side perturbatively. In this analysis, we assume 
that the external momenta are close to their mass shell, — ml ~ mi,A and py ~ A^, and 
that the momentum transfers scale as pb ■ Pv ~ P_b ■ P 7 ^ Pv ' P'r ~ where A remains 
hxed in the heavy-quark limit. The correlator is then expanded about the heavy-quark limit. 
Perturbative factorization relies on the assumption that if one hnds that the double spectral 
density of the correlator in (5) with respect to the variables p^ and py has certain factorization 
properties to a given order in A/mi, and to all orders in perturbation theory, then the same is 
true of the amplitude on the right-hand side of the reduction formula. For operators such as 
Q\ 2 containing charm quarks, we need to specify how the charm-quark mass is treated in the 
heavy-quark limit. We take the limit holding the ratio me/mi, fixed. 


2.1 Diagrammatic analysis and momentum regions 

The correlator in (5) is a function of Lorentz-invariant scalar products of the external momenta. 
However, to obtain its expansion in powers of l/rrif, it is advantageous to introduce reference 
vectors. We introduce a unit four-vector in the direction of the B meson and a light-like 
vector in the direction of the outgoing vector meson, and dehne 

fi. = J_ h.-- itiP (6) 

n ■ V \ n ■ V J 

so that = 0 and n-n = 2. We decompose all momenta into their light-cone components, 

p^ = n ■ p — + n ■ p — + p^^ = pI + p'/L + p'/^. (7) 

Note that our dehnition implies v±_ = 0. Working in dimensional regularization, we then 
employ the strategy of regions [26, 27] to expand diagrams about the heavy-quark limit. 
With this technique, the integrands are expanded in a number of different momentum regions. 
The expansion of the full integral is recovered after integrating the expanded integrands and 
summing the contributions from the different regions. 

The momenta in the relevant regions differ by the scaling of their components (p_|_,p_,px). 
Not surprisingly, momentum regions in which a loop momentum has the same scaling as an 
external momentum give a non-zero contribution. These regions are: 

soft: (A, A, A) mi, , 

n-collinear: (A^,l,A)mf,, 

where we introduce a dimensionless expansion parameter A ~ A/mi,- The loop momentum 
scales like ps — mi,v in the soft region, and in the same way as pv in the n-collinear region. 
An n-collinear region does not appear, since p^ = 0 (it would be present for B if 

^ addition to these scalings, regions with p^ ^ A^ arise: 
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hard: (1,1,1) rrib , 

n-hard-collinear: (A, 1, V^) mi, , 
h-hard-collinear: (1, A, V^) mi, . 

The presence of two large perturbative scales — the hard scale ~ ml and the hard-collinear 
scale ~ mfeA — manifests itself in the factorization formula: the hard-scattering kernels 
have the schematic form T ~ C* ® J, where C and J include hard and hard-collinear contri¬ 
butions, respectively, and the symbol denotes a convolution over momentum fractions. 
The n-hard-collinear region arises in diagrams where the photon attaches to the spectator. 
Offshell propagators also appear; for example, a momentum scaling (A, l,A)mfe arises when 
a soft and n-collinear momentum flow into the same vertex. This will be discussed in more 
detail in Section 2.3. Finally, the following low-energy region appears in the expansion of the 
diagrams [25]: 


n-soft-collinear: (A^, A, A^/^) m;, . 

This momentum region appears in interactions with soft and collinear lines. The scaling of its 
components is the largest compatible with both the soft and the collinear scaling: {ps+Ps<^‘^ ~ 
pI and (pc + Psc)^ ~ Pc- Note that {ps +PcY ~ pL so that collinear lines cannot emit or absorb 
soft momenta and remain collinear, or vice versa. Since it is the only low-energy interaction 
connecting the soft and collinear sectors, proving factorization to a given order in A amounts 
to showing that there is no contribution from the soft-collinear region to this order. 

Figure 1 shows three typical contributions to the decomposition of the correlator (5), for 
the chromomagnetic operator Q^. The three rows in the hgure illustrate the soft-overlap 
(A), hard-scattering {B), and spectator-emission {C) mechanisms. The two-step matching 
procedure QCD —SCETj —s- SCETn is described in the following Sections 2.2 and 2.3. 

Note that the soft-collinear region has p^ ~ /mi, -C A?. It has sometimes been argued 
that it is “unphysical” to allow for momentum regions with p^ parametrically below A^, since 
non-perturbative effects would modify physics below this scale, and that it would therefore be 
more natural to perform the perturbative factorization analysis with a hard infrared cut-off 
in QCD. Since the key point in factorization proofs is precisely to show that such infrared 
regions are either absent or cancel in the sum over diagrams, simply ignoring such modes 
is clearly not an option. If one chooses to introduce an infrared cut-off in QCD, the proof 
of factorization becomes equivalent to the demonstration of insensitivity to this regulator. 
However, it is difficult to introduce such a cut-off in a gauge invariant way,^ and it is also 
doubtful whether the diagrammatic analysis with a cut-off can be reformulated in effective- 
theory language. Since the messenger helds do not contribute (by dehnition) to factorizable 
quantities, and since non-factorizable quantities are categorized as non-perturbative, nothing 
is gained by removing these helds in favor of an infrared cut-off. One of the advantages of the 
effective-theory approach in dimensional regularization is precisely that the analysis can be 
performed without explicit momentum cut-offs. 

^The only known way is to quantize in a finite volume and use twisted boundary conditions for the gauge 
fields to eliminate the zero mode. 
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QCD 



SCETi 



SCETii 



Figure 1: Three QCD Feynman diagrams for the contributions of Qs and their leading-order 
representation in the effective theory. The double line denotes the heavy-quark held. The 
dashed lines denote hard-collinear helds in SCETj and collinear helds in SCETn. Solid lines 
in the effective-theory diagrams denote soft helds and the dotted line denotes a soft-collinear 
held. 

While it is easy to see that all of the above regions are required to obtain the expansion of 
the correlator diagrams, we do not have a proof that they are sufficient.^ Two-loop applications 
in similar kinematic situations [28] suggest that no additional regions are needed. The above 
list of momentum scalings is natural in that it contains all onshell modes whose components 
n-p and n-p scale with powers of A equal to the scaling of the components of external momenta. 

Finally, let us note that the analysis of regions presented above assumes exactly massless 
light quarks. A systematic inclusion of quark mass terms presents a challenge, since the mode 
structure in the low-energy theory is then drastically altered. For instance, including 0{A.) 
masses would eliminate the soft-collinear mode, but the resulting diagrams for the soft and 
collinear regions would no longer be separately well-dehned in dimensional regularization, 
requiring additional unconventional (e.g., analytic) regulators. We will return to this issue in 
Section 4.3 and address the more modest question of the leading corrections for light-quark 
masses niq A. We argue that contributions linear in the light mass may be absorbed into 
the hadronic parameters appearing in the factorization formula, while any terms that could 
potentially spoil factorization appear Erst at quadratic order. 

^The same is true for traditional diagrammatic factorization proofs. Additional momentum regions could 
invalidate the analysis also in these cases. 
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2.2 Intermediate effective theory: SCETj 


In the construction of the SCET Lagrangian, an effective-theory held is introduced for each 
momentum region. The integrands of the QCD Feynman diagrams expanded in the various 
regions are then reinterpreted as arising from the Feynman rules of the effective theory. Fur¬ 
thermore, in order to ensure that the amplitudes are appropriately expanded in momentum 
space, a derivative (“multipole”) expansion is performed in the effective action. 

Note that if we had chosen the momentum of the vector meson to be n-hard-collinear, 
then only the hard, n-hard-collinear, h-hard-collinear, and the soft region would appear in 
the expansion of (5). It is simpler to hrst consider the situation where we count the external 
momenta in this way and to introduce helds only for these regions. This is illustrated by the 
middle column of Figure 1. The corresponding effective theory, called SCETi, describes QCD 
at or below the hard-collinear scale, and contains the quark and gluon helds 

~ ~ (A, 1, v^), ^71^ ~A^/^ A^~(l,A,yA), 

( 8 ) 

qs ~ A^/^, ~ (A, A, A), h ~ 

The hard-collinear quarks are described by two-component spinors satisfying ~ 0- 

We have indicated in (8) the scaling of the held components, which can be derived from 
the scaling of the corresponding propagators. No held is introduced for the hard region, as 
this contribution will be absorbed into the Wilson coefficients of operators in the ehective 
theory. As was shown in [5], for diagrams involving only a single type of hard-collinear held, 
in dimensional regularization there is no hard contribution in the pure QCD sector, since 
the corresponding diagrams are scaleless and vanish. The ehective Lagrangian can then be 
constructed exactly, to all orders in perturbation theory. This was done to next-to-next-to- 
leading power in [29]. The same may be done in our case with two types of hard-collinear 
helds, which we denote generically as (phc and 0^^. In fact, the Lagrangian for this case is 
simply 

^ ( 9 ) 

where Cg is given by the HQET Lagrangian for heavy quarks, and by the restriction of the 
QCD Lagrangian to soft momentum modes for light quarks and gluons. Che denotes the 
remainder containing the hard-collinear helds and their interactions with the soft helds, and 
Cj^ is obtained from Che by interchanging n and n. 

Note that we did not write down a Lagrangian containing interactions with both n- and n- 
hard-collinear helds. By the Coleman-Norton theorem [30], pinch singularities can only occur 
in momentum conhgurations that can be interpreted as classical scattering processes. By 
momentum conservation, this cannot occur in interactions with both n-hard-collinear and n- 
hard-collinear particles unless both types of particles are present in the initial and hnal states. 
An example is illustrated in Figure 2, where the momentum k is restricted to the region k ~ 
(A, A, V^). When this skeleton diagram is inserted into loop diagrams, an integration over k in 
this region involves denominators of the form {phc±+k±)'^+n-phen-k and {pj^j_+kj_)‘^+n-p^n-k. 
The contour integrals in n ■ k {n ■ k) vanish if the external particles all have the same sign 
of n ■ Phc {n ■ Pf^), and such exceptional conhgurations are therefore not relevant in cases 
where collinear particles are present only in the hnal state. The absence of such exceptional 




Figure 2: Example of an exceptional momentum configuration giving rise to an interaction in 
SCETi involving he and he fields. 

momentum configurations is encoded automatically in the usual strategy of regions applied to 
the perturbative expansion of the correlator (5) in dimensional regularization. For an iV-loop 
diagram, this strategy assigns onshell momentum scaling to N internal lines, with the scaling 
of all remaining lines hxed by momentum conservation. The full amplitude is recovered by 
performing this assignment in all possible ways, with an unrestricted integration over the N 
onshell loop momenta. With these rules, an isolated offshell line such as in Figure 2 (with 
momentum phe + + k) cannot occur. 

It is convenient to use the SCETj operators to classify the different mechanisms through 
which the decay B ^ V'y can proceed. In contrast to the Lagrangian interactions, there 
are hard matching corrections to the weak-interaction operators in the effective theory. After 
performing the matching of QCD onto SCETj, the remaining problem is to examine in each 
case all possible SCETn operators that can result. We now turn to this problem and discuss 
the issues involved in integrating out the hard-collinear components of the SCETj helds. 


2.3 Final effective theory: SCETn 

Counting the external momenta as collinear, instead of hard-collinear, the full list of regions 
in Section 2.1 needs to be considered. The dynamical helds in this case are the collinear and 
soft helds 

C~A, A^~(A^,1,A), ~ A^/^, A^~(A, A,A), h ~ A^/^, (10) 

as well as the soft-collinear quark and gluon helds 

yl" ~(A^A,A=>-'2). (11) 


The collinear and soft-collinear helds are again described by two-component spinors satisfying 
ii^c = ‘i'Q = ^- The small-component projection of the soft-collinear fermion held, satisfying 


= 0, is given by 


a = — 


fji 1 
2 in ■ Dsc 


'P..J ■ 


( 12 ) 


In SCETn, both the hard and hard-collinear contributions are absorbed into the Wilson coef- 
heients of the operators built from the above helds. The hard-collinear contributions appear 
in the matching step from SCETj onto SCETn. As in the hrst matching step, the pure QCD 
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(A, A, A) 


(A^1,A) 



(A,1,A) 


Figure 3: Example of a SCETj diagram relevant for B Vj decay with an offshell-collinear 
gluon. The same contribution is also depicted in the second line of Figure 1. 


part of the effective Lagrangian can be obtained exactly. It was constructed at next-to-leading 
power in [25]. 

In preparation for the discussion of general operator bases to be considered in Section 3, 
we review here the procedure employed in integrating out the helds at the hard-collinear scale. 
We begin by restricting attention to the sector of the SCETi Lagrangian (9) involving n-hard- 
collinear and soft helds. Just as matching QCD to SCETi involved decomposing the QCD 
helds into their hard, hard-collinear, and soft components, so matching onto SCETn involves 
the decomposition of SCETi helds. In particular, a generic hard-collinear held is decomposed 
as 

4>hc 4>hc + 0c + 0 OC •} (13) 

where the hard-collinear, collinear and ohshell-collinear momenta scale as (A, 1, V^), (A^, 1, A), 
and (A, 1, A), respectively. The latter momentum scaling arises from the combination of onshell 
soft and collinear momenta, see Figure 3. The hard-collinear and ohshell-collinear helds are 
integrated out in passing from SCETi to SCETn, just as the hard components of the QCD 
helds were integrated out in passing from QCD to SCETi. Similarly, a generic soft held is 
decomposed as 

0s ^ 0s T 0SC ) (1^) 

where the soft and soft-collinear momenta scale as (A, A, A) and (A^,A,A^/^). In contrast to 
SCETi, where the soft helds are dehned to contain all modes below the soft scale, the soft helds 
in SCETn are dehned to contain strictly soft, and not soft-collinear modes. This interpretation 
is mandated by the appearance of a new region in SCETn. If we were to work with explicit 
cut-ohs, the soft helds would be required to have n ■ ps of order A (thus excluding modes 
with n ■ ps oi order A^), and collinear helds would be required to have n ■ pc oi order unity 
(excluding modes with n-pcoi order A). The situation is analogous to the passage from QCD 
to SCETi, where the hard-collinear region is dehned to contain strictly hard-collinear, and not 
soft, modes. 

We now expand the SCETi Lagrangian (9) using the decompositions (13) and (14). We 
split the Lagrangian into two parts, one containing the “light” degrees of freedom present in 
the low-energy theory, and one containing “heavy” modes that are to be integrated out. In 
the hrst part, Dught, we collect all terms that contain only helds that are part of SCETn: soft. 
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collinear, and soft-collinear fields. Ziught has been derived in [25] and is required through 0{\): 

r — r r r -u -i- -u 

-‘^hght \ i—c T )—sc 1" -‘-s+sc ' ’—c+sc “T • • • 

_ £{0) £(1/2) £{1) /i['\ 

Terms with both soft and collinear fields appear at subleading power in the decomposition of 
the SCETi Lagrangian, both directly in £iight, and via induced interactions after integrating 
out offshell modes in £heavy below, as discussed in [6]. However, such interactions are not 
relevant to our analysis, as they do not appear in the expansion of the correlator (5). More 
generally, they are absent in cases where collinear particles are present only in the final state, 
by the same reasoning as for the terms with both n-hard-collinear and n-hard-collinear fields 
in the decomposition of the QCD Lagrangian in Section 2.2. 

In the remaining part of the Lagrangian, /iheavy, we collect all terms that involve at least 
one hard-collinear or offshell-collinear field, which will be integrated out in the construction of 
SCETn. For simplicity, in the discussion of Lagrangian terms involving such “heavy” modes, 
we work in light-cone gauge n ■ H = 0 for the fields descending from A^c (he., hard-collinear, 
collinear, offshell-collinear), and n - H = 0 for the fields descending from Ag (i.e., soft and soft- 
collinear). To fully separate the different scales, interactions involving fields with different 
momentum scaling must be multipole-expanded and the offshell-collinear fields ^oc and A^c 
integrated out. The remaining onshell fields can be assigned a definite power counting, and 
the offshell fields are expressed in terms of a series (ordered in A) giving the possible branchings 
into these onshell fields. For interactions of collinear with hard-collinear fields we have 

(t)hc{x)(l)c{x) = (l)hc{x) 0c(a:+) + a;([(9^0c(a:+) + {x^-d^ + ]^x^^x''^d^dy)(t)c{x+) + ... . (16) 

Similarly, for soft and hard-collinear fields ((>hc{x)(j)s{x) ~ (l)hc{x)(l)s{x-), while for soft and 
collinear fields, (j)s{x)(j)c{x) ~ (l>s{x- + x±)(f)c{x+ 4- xj,). 

We first expand £heavy in powers of A. We begin with the tree level case (i.e., neglecting 
interactions involving onshell hard-collinear fields), where we will find that the solutions for 
the offshell fields scale as ^oc A^/^, and n ■ Aoc ~ A^. We will then consider the 

inclusion of hard-collinear fields, finding that the solutions for the offshell fields in this case 
start at one power lower in A: 

C _ ty(l/2) , f(l)\ , t(3/2) , 

S,oc — \^S,oc ‘ Soc J ‘ Soc . . . , 

■4«x=(Lk>+4»i)+42f + .... (17) 

n-A„^{n- ylW + n ■ + n • + ... . 

The terms in parentheses only appear when branchings into hard-collinear fields are included. 

The tree-level Lagrangian begins at 0{X), and for a complete matching at leading power 
we require terms through (9(A^). With the inclusion of hard-collinear fields, the Lagrangian 
begins at one power lower in A: 

r _ { /’(O) I ri^lX) 5 ^(1) I 

-^heavy — I -^heavy “r '^heavy J ' '^heavy “r . . . . 
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Omitting the terms involving hard-collinear fields for the moment, the leading fermion La- 
grangian reads 

^htLy = ■ <9+ {loJ^^g4^cLqs + + h.c.} , (19) 

and solving the equation of motion yields 

Similarly, from the gluon terms in we find 

= - /■ , {r9.7^r;, + h.c.}. 

in - dm - o ^ 

= ( 21 ) 

Having found the leading terms, these solutions may be substituted back into £heavy, and the 
process iterated at the next power in A. The complete list of SCETn operators at leading 
power requires also and these terms themselves involve and n ■ A^oJ^\ The 

tree-level expressions were obtained in [11].^ 

Beyond tree level, we must consider the branching of offshell-collinear fields into two or 
more onshell hard-collinear modes. At each order in A, we expand in powers of the coupling 
constant g. Only factors of g associated with hard-scale (i.e., hard-collinear and offshell- 
collinear) gluons are included in this expansion. Anticipating that and ~ 

g, contributions to /Iheavy involving offshell-collinear fields begin at 0{g‘^). For the fermion 
Lagrangian, 

^heLy = ^hc + (^gn- Ahc + g4hc±:r^i^±'^ ihc + h.c.j 

+lhc^ n ■ AW T^g4oc±^ + g4oJ± ihc + 0{g^). (22) 

Solving the equation of motion yields 

(n • xlxo + ^0 + 0(jt. (23) 

in ■ o \ in - o J 

From the gluon terms in we find in the same manner 


in - din - a y 2 \ in - o in - d J 

^Expressions for our ^(T^t(2),(5/2) given by in [11]. In the decomposition of SCETi fields at 

tree level in this reference, “Ac” refers to what we call “oc”. Explicit expressions are given there for the slightly 
different quantities ■j/)(3)>(4),(5)^ which include contributions from the soft field qs and from the small-component 
hard-collinear field ghc in SCETj. These terms are not part of ^oc) in particular, the terms containing from 
(m • 5)“^(z.^c^-I-(/4s-L)(fl/2)^c from and (m • 5)“^g4ic± Hf^/2)Cc from should not be included 
in ^oc- 
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+ \[‘Kc^,in ■ dn ■ Akc] - [n in ■ 9/1^1 “ [■4(ic±^,ia'I^L “ j + 0(9^) , 

n ■ + r“{/.o|r“&, j + C>(/), (24) 

Note that the hard-collinear helds appearing on the right-hand side in (23) and (24) must 
have total transverse momentum of order A, even though the individual helds have transverse 
momentum of order A^/^. This constraint is automatically enforced by the scaling of external 
soft and collinear momenta in the evaluation of diagrams corresponding to these interactions. 
Having derived , and n ■ to the desired order in g, these solutions can be 

substituted back into £heavy, and the procedure iterated at the next order in A. This process 
can be carried out to any order in the power expansion. For offshell-collinear helds branching 
into hard-collinear helds, the complete list of SCETn operators at leading power requires 

fl) (3/2) fl) f3/2') 

also ^oc and ^oc , which themselves involve A)^^^ and n ■ Aoc ■ Since further subleading 
Lagrangian interactions are required to convert the remaining onshell hard-collinear helds 
into soft and collinear partons, such operators are required only up to one power in A lower 
than in the tree-level case. 

Substituting the expressions (17) for the ohshell helds into £heavy in (18), and inserting 
appropriate gauge strings to relate the expressions in light-cone gauge to those valid in an 
arbitrary gauge, yields the hnal result for the decomposition of the SCETj Lagrangian in 
the n-hard-collinear sector. The sector of SCETj involving h-hard-collinear modes may be 
treated similarly. It is convenient to treat the photon as being an n-hard-collinear held in 
the intermediate ehective theory and as an n-collinear held, ~ A, in SCETn. The hnal 

results are independent of any power counting assigned to this held, since we work to hrst 
order in the electromagnetic coupling. As discussed in Section 2.1, no other n-collinear helds 
appear in the B ^ V'y analysis. The decomposition of the SCETj Lagrangian in this sector 
is obtained from (18) simply by replacing n h and dropping all other n-collinear helds. 
The same manipulations as in the previous case yield the hnal result for the explicitly gauge- 
invariant and multipole-expanded Lagrangian in the sector involving n-hard-collinear and soft 
helds. The matching of SCETi onto SCETn is completed by substituting the solutions for 
the ohshell helds into external current operators, and integrating out the remaining onshell 
hard-collinear helds. The hard-collinear modes do not contribute additional renormalizations 
to the relevant part of the low-energy QCD Lagrangian, but result in non-trivial matching 
conditions for external currents. This matching is discussed in detail in Section 3. 

As an illustration of the passage from SCETj to SCETn, we may consider the represen¬ 
tation of operators contributing to form-factor matrix elements. The leading-power SCETi 
current operators are of the schematic form ^hch- Using the decomposition (13) and enforc¬ 
ing momentum conservation to drop the term involving a single hard-collinear (he) held, the 
mapping onto SCETn operators is given by: 

^hch ^ch + ^och . (25) 

Expanding the solution of the equation of motion for the held ^oc as in (17), we then have 
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diagrammatically at tree-level; 



s c s s s c 



s c s c s c 


( 26 ) 

The contribution of the leading operator ^ch ~ receives an additional suppression when 
inserted into correlator diagrams (analogous to Figure 1) for the form factor, because these 
diagrams will always involve soft-collinear quark lines. For example 


A5/2 C A5/2 C 



Ai/2 sc A A sc ^1/2 


With the additional A^/^ suppression from subleading Lagrangians terms in £iight describing 
the coupling to soft-collinear quarks, the contribution ends up being of order A"^, the same 
order as the contribution of Another example of this additional suppression relating 

to i? —1^7 is illustrated by the first row in Figure 1, where the operator appears in 
combination with subleading interpolating current operators for the initial- and hnal-state 
mesons that contain soft-collinear quarks. We will discuss this in more detail in Section 3. 

The operators from ~ A^ and their leading-order matching coefficients were given in 

[12]. Additional terms arise at leading power from for flavor-singlet hnal-state mesons 

and have not been shown in (26). As discussed in more detail in Section 3, these terms 
contain collinear gluon helds in place of the collinear fermion bilinear. For the remaining 
terms, A^ gives rise to soft-overlap contributions in the havor-singlet case, connected 

with terms arising from Similar to the operators in ~ A^/^ receive an 

endpoint suppression, and as a result do not contribute at leading power. 

The subleading SCETj operator ^hcd±h may be decomposed in a similar way: 

Chcd±h ^ ^cd±h + ^ocd±h . (28) 

However, in each case the d± derivative gives an additional 0{\) suppression relative to 
the operators in (25). The remaining form-factor contributions of leading power arise from 
subleading SCETj operators of the form ^hcAhc±h. The decomposition in this case is 

^hc-^hclh ^ ^hc-^hclh ~\~ ^c-^oclh ~\~ ^oc-^cJ-h ^oc-^ocA-h . (^ 9 ) 
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Again, from the expansions (17) we find at tree-level; 



c 

(gmAZ'f’h = OiXO'^) (30) 

The operator ^ contributes hard-scattering contributions of leading power, while 

Aci_h ~ A*^ contributes only for flavor-singlet hnal-state mesons. Because of an additional 
A^/^ endpoint suppression, icAci.h ~ A’^^^ cannot contribute at leading power. Likewise, the 
contributions of h A®/^ are also power-suppressed. 

This procedure can be extended beyond tree-level by integrating out the (onshell) hard- 
collinear modes. For instance, for the first term on the right-hand side of (29), 




he AF2 he aF2 



There are also contributions where offshell-collinear helds branch into hard-collinear modes. 
For instance, from the second term on the right-hand side of (25), 
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In (32) we have displayed a contribution involving the solution for the offshell-collinear gluon 
field substituted back into /Iheavy (22). While straightforward in principle, these 

examples illustrate the nontrivial nature of the SCETj to SCETn matching. Instead of ex¬ 
plicitly integrating out the hard-collinear modes, in Section 3 we will arrive at the complete 
SCETn operator basis using only general properties of the decomposition of SCETj operators. 
The general form is required both for explicit computations, and to demonstrate factorization 
properties, such as the decoupling of leading-power soft-collinear interactions. In contrast to 
SCETi, the power counting for operators in SCETn cannot be deduced simply by inspection 
of the field content. This is illustrated by the third line of (26): the two operators with an 
additional gluon held turn out to be of the same order as the four-quark operator, due to the 
non-localities introduced by integrating out the hard-collinear modes, of virtuality ~ rribA. 
These non-localities manifest themselves as inverse partial derivatives, counting like [6]. 
The appearance of these derivatives is manifest in (20) and (21), the solution of the equations 
of motion for the offshell-collinear helds. In order to proceed, we require a set of rules that 
can restrict the appearance of such factors and, more generally, allows us to write down the 
most general SCETn operators. 

3 SCET representation of the weak Hamiltonian 

As discussed in the previous section, the QCD part of the low-energy effective theory does 
not receive matching corrections and can be constructed exactly. This is not true for the 
operators in the weak Hamiltonian. For this case we proceed in the usual way: we write down 
all operators with the correct quantum numbers built from the available fields, and perform 
perturbative matching to the desired order. Our goal in this section is two-fold: to hud the 
SCETn operators that contribute to B Ey decay at leading power, and to construct the 
SCETi operators which match onto these SCETn operators. The utility in identifying the 
SCETi operators lies in the fact that the soft-overlap contributions can be isolated already at 
this stage, before further decomposition into SCETn fields. The two-step matching procedure 
is also required for the resummation of large perturbative logarithms, which we address in 
Section 5. Using building blocks defined below, it is straightforward to write down all SCETi 
operators that can contribute up to a given order in A. The situation is more complicated 
for the SCETn operators: when integrating out hard-collinear modes, inverse derivatives n- d 
on the soft helds appear [6], counting as A“^. Despite the presence of such derivatives, we 
will see that only a hnite number of operators can appear to a given order in A. A second, 
practical difficulty is that the leading SCETn operators are of a much higher order in A than 
the leading SCETi operators. For example, treating the photon held as a hard-collinear held 
in SCETi and as a collinear held in SCETn, the leading operator in the intermediate theory 
contributing to B ^ V'j counts as A^/^, while the leading operators in the hnal ehective theory 
count as A^. Power counting alone does not strongly constrain the possible SCETi operators, 
and would leave us with a very large number of operators in the intermediate theory, most of 
which would turn out to be irrelevant upon matching onto SCETn. We will hnd that counting 
the mass dimension of the SCETi operators leads to much stronger restrictions. 
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3.1 Building blocks 

A characteristic feature of SCET is that derivatives of the (hard-)colhnear helds corresponding 
to large momentum components are unsuppressed, and operators with an arbitrary number of 
such derivatives can appear at the same order in the power counting. To account for this, the 
operators are allowed to be non-local along a light ray: for example, the SCETj representation 
of a QCD operator at position x can contain the hard-collinear helds (phdx + sn), 0^(a;-l-rn). 
The Wilson coefficients of the operators are then functions of the light-ray variables (r and s 
in our example). 

To obtain gauge-invariant operators, the helds at diherent points on the light-ray must be 
connected by light-like Wilson lines (deviations of such Wilson lines from the light cone can be 
expanded and appear as power-suppressed operators). Instead of inserting these Wilson lines 
for each operator, it is simpler to work with building blocks [6, 31] obtained by multiplying the 
helds by Wilson lines which run along the light-ray to inhnity. These building blocks will be 
invariant under hard-collinear gauge transformations in SCETj, and under soft and collinear 
gauge transformations in SCETn. We choose to work with building blocks that have simple 
transformations, but are not invariant, under soft and soft-collinear gauge transformations 
in SCETi and SCETn, respectively. Purely gauge-invariant quantities may be obtained by 
introducing additional soft or soft-collinear Wilson lines, but this will not be necessary for our 
arguments, and would require the appearance of residual Wilson-line factors in SCET current 
operators. The building blocks dehned here are also easier to work with when performing 
explicit loop calculations. Thus, for SCETi, we introduce the helds 

‘Xhcix) = Wl{x)^Ux), (33) 

KSA = ^Viix) [WlXx) WUx)] + y [Wl(x)gn ■ A.{x.)W,Ux) - gn ■ A.ix.)] , 
with = idfg -|- gAhc and Wilson line 

Whc{,x) = P exp i^g j dsn- Ahc{x -t- sn)^ . (34) 

Note that n ■ A'^^{x) = 0. The building blocks for the hard-collinear helds in the opposite 
direction, ^j^{x) and Aj^{x), are obtained by interchanging n and n (and X- —> a;+) in the 
above expressions. 

The building blocks of SCETn are dehned in an analogous way. In this case the role of 
the soft helds is played by the soft-collinear helds, and both the soft and the collinear helds 
are supplied with Wilson lines: 

Xc{x) = W^{x)^c{x ), 

77 ^ 

A^^{x) = Wl{x) [iD^{x) Wc{x)\ + — \wl{x)gn ■ Ascix-)Wc{x) - gn ■ A^dx-)] , 

(35) 

Qs(a;) = S'](a;)g^(a;), XCs(x) = Sl(x)h(x) , 

A^(x) = S^x) [iB^(x) S' 5 (a;)] -F ^ [Sl{x)gn ■ Asc{x+)Ss{x) - gn ■ Asc{x+)] . 
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The collinear Wilson line Wc{x) is defined in the same way as Whc{x) in (34), except that it is 
constructed with the collinear instead of the hard-collinear gluon field. The soft Wilson line 
is ^ 

Ss{x) = P exp (^g j dtn ■ As{x + tn)^ . (36) 

For a detailed discussion of the gauge transformation properties of the SCETn fields and the 
construction of gauge-invariant building blocks, we refer the reader to [25]. Similar to the 
building blocks for the n-hard-collinear fields, Xj^{x) and Aj^{x), we will also need SCETn 
building blocks for which n- and h-directions are interchanged. The only collinear field in the 
h-direction is the photon field. However, we will need the associated soft building blocks with 
Wilson lines in the h-direction and will denote them by Qs(a;), “Ksix) and As{x)^ 

Arbitrary SCET operators are obtained by combining the above building blocks. In prod¬ 
ucts involving different momentum modes, a derivative expansion of the fields has to be per¬ 
formed [5, 25]. The expansion for fields in SCETj is as in (16), and for SCETn we have 


(l)s{x)(l)c{x) (t>s{x_ + X^)(l)c{x+ + X±) 

+ x+ ■ ds(t)s{x_ + Xi_)(t)c{x+ + x±) •+• (l)s{x_ + x±)x_ ■ (9c0c(a;+ -f- x^) + ... , (37) 

and similarly 0s(a;)05c(a:) ~ 0s(a;)05c(a:+), 0c(a^)0sc(a^) ~ 4>c{x)(j)sc{x_). Eor our leading-power 
analysis the derivative terms can be dropped, and we will suppress the x-dependence of the 
various fields in the following. 

3.2 Operators in SCET 

We now present a general procedure for matching generic SCETi operators onto SCETn. The 
SCETi operators are products of soft fields and hard-collinear fields in the n- and h-directions; 
schematically we may write 

O = [“soft”] X [“n-hard-collinear”] x [“h-hard-collinear”]. (38) 

Because the SCETi Lagrangian (9) decomposes into the two hard-collinear sectors, each of 
the three brackets can be treated separately and they match as follows: 

[“n-hard-collinear”] —> [“n-collinear”] x [“soft”] x [“soft-collinear”], 
[“h-hard-collinear”] —x [“soft”] x [“soft-collinear”], (39) 

[“soft”] —^ [“soft”] X [“soft-collinear”]. 

Physically, the reason that the sectors match separately can be understood by picturing the 
decay process: at a certain time, the heavy-quark decays into two energetic partons flying 
in opposite directions. Each of these two particles can subsequently emit soft and collinear 

^As{x) contains a soft-collinear gluon field in the opposite direction. However, since there are no collinear 
quark or gluon fields in the h-direction, and since the messenger fields only contribute in exchanges between 
soft and collinear particles, this region does not contribute in H —> Hy. 
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particles, but the energetic particles from opposite directions cannot annihilate each other. 
This physical picture is formalized by the Coleman-Norton theorem. The fact that the soft 
sector matches separately follows because the soft fields are not integrated out in the transition 
to SCETii. Each soft held in SCETi is simply replaced by the sum of a soft and a soft-collinear 
held in SCETn, see (14). 

One complication is that the individual sectors are generally not invariant under soft gauge 
transformations, while our SCETn building blocks are invariant. In most cases we can avoid 
matching non-invariant operators by gronping the “soft” bracket together with either the “n- 
hard-collinear” or “n-hard-collinear” bracket. In the general case, we can introdnce soft Wilson 
lines to make each sector gange invariant and remove them after the matching is completed. 
We shall come back to this point in Section 3.2.2. 

3.2.1 Current operators 

We hrst discuss the simplest case, namely SCETi operators of the form 

O = [“soft”] X [“n-hard-collinear”] x , (40) 

where we use the schematic notation of (38) for the special case where only the photon field 
appears in the “n-hard-collinear” bracket. These operators arise when the photon is emitted 
from one of the cnrrent qnarks, bnt onr discussion does not depend on this fact. The analysis 
for this case is identical to that for the cnrrent operators dehning heavy-to-light form factors. 
The constrnction of the general SCETn operator basis relevant at leading power has been 
performed in [11, 12]. We now rederive these resnlts as a preparation for the general case, 
and to introduce our method. We start by writing ont a list of the lowest-dimension current 
operators in SCETj. By momentnm conservation, the operators must contain at least one 
hard-collinear held. We will see below that it is most convenient to classify SCETj operators 
according to mass dimension, rather than power connting in A. Up to dimension hve, we End: 

d = 3:Xn,rh (J^), 
d = 4 : Xhc Ahci. T' h (J^), Qs A.hc± T'^ h , 

Xhc dhc± h , Xhc Ds± T' h , 
d 5 . Xhc Ahcl. dI/ic_L T h , Ahcl. dI/ic_L E Yi , 

Xhcn-dhcU- Ahc r' h , Xhc Xhc Xhc r' h , (41) 

n-dhc2 

Xhc _ Xhc Qs r" h , qsfi- dhc n ■ Ahc E" h , 

n- dhc 2 

Qs ^ ' dhc ^/ic_L /i , .... 

The symbols in parentheses, and J^, anticipate the notation to be introdnced in Section 4 
for the relevant SCETi operators. We do not display transverse Lorentz indices or color indices; 
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the former may be contracted with the metric and epsilon tensor in the transverse plane, 

gT = • (42) 

We use the convention = —1. In the schematic notation of (41), it is understood that 
hard-collinear derivatives can act on any of the hard-collinear helds in the operators, and 
similarly for soft derivatives. We do indicate the Dirac structures that can occur in the above 
expressions, using the following Dirac matrices that are invariant under the “boost” > an^ 
and a~^n^\ 

r' = {1,75,7!}, 

r" = r'u{#^,7!75,7!7l-7!7!}. 

The sixteen matrices ifiT', fiV, and T" form a Dirac basis. We only consider boost-invariant 
operators; such a choice is always possible and is also natural because the operators we repro¬ 
duce with the effective theory are independent of the reference vectors and Operators 
that are not boost-invariant can be eliminated in favor of invariant operators obtained by 
multiplying them with an appropriate number of derivatives n ■ dhc- The last hve operators 
of dimension hve are examples where such derivatives have been included. The presence of 
these derivatives can be compensated by the Wilson coefficients of the non-local operators, for 
example 



j ds C(s) n ■ dhc4>hc{x + sn)(j)s{x) 



4>hc{x + sn)(l)s{x). 


(44) 


We did not allow for operators which explicitly involve the vector in (41) because it can 
be eliminated in favor of n^, cf. (6). Furthermore, we have used the projection properties of 
the spinors Tj>h = h, i/i'Xhc = 0 to eliminate occurrences of ^ h. For instance, the third operator 
for d = 5 is obtained by rearranging an operator with d = 4: 

n/i 

2n ■ V 

^hc ' dhc ^ ' d\.}ic F h . 


n ■ V 


■ n ■ Ahr h = 


Xhr F'n ■ Ahr h 


n ■ V 


(45) 


Ln.A^h = X^T’ 


In the second line, we have absorbed a factor {in ■ dn ■ v)~^ into the Wilson coefficient of the 
operators, as in (44). To minimize the list of possible SCETj operators appearing with a given 
dimension, it is convenient to always make use of such rearrangements.® 

®There are additional constraints arising from the independence from the reference vectors. Requiring 
complete reparameterization invariance, also under ^ n^ + e([, yields relations linking operators of different 
orders in the power counting. Since we are concerned only with the leading order, such transformations will 
not be relevant to the present discussion. 

^Another possibility would be to group the appearing on the left-hand side of (45) together with the 
heavy-quark field. Since the heavy-quark does not participate in the matching of SCETi onto SCETn, the 
general SCETn operator is given by examining the matching of a boost non-invariant operator containing 
n • Ahc- This approach is taken in [11], and using a larger set of building blocks, such operators can be shown 
not to contribute at leading power. 
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Table 1: Boost-invariant building blocks for SCETn operators, with their dimension d and 
order [A] in the power expansion. Soft derivatives dg can act on any soft held in the operator, 
collinear derivatives dc on any collinear held. T' and T" are dehned in (43), while and e'^ 
are dehned in (42). Additional building blocks are obtained by hermitian conjugation or by 
replacing Qs with "Kg. 


For contributions arising at leading power, the SCETi operators should not contain soft 
helds in addition to those found in the hnal SCETn operators. Such soft helds would not 
participate in the matching, and result in a power suppression relative to the corresponding 
operators without the additional soft helds. Similar arguments apply to power-suppressed soft 
or collinear derivatives; an explicit example of this ehect was mentioned in (28) of Section 2.3. 
Thus only the hrst two operators for d = 4 can be relevant for our leading-power analysis. We 
have also not listed operators of any dimension containing additional factors of n ■ Dg/n ■ dhc- 
The ellipsis for d = 5 denotes similarly irrelevant terms. 

In order to construct all SCETn operators up to a given power, we work with the set 
of building blocks in Tables 1 and 2, which are invariant under soft and collinear gauge 
transformations. Again, we choose to work with boost-invariant quantities. We will begin 
by using Table 1 to describe the matching onto SCETn operators corresponding to “typical” 
momentum conhgurations in which the partons in the initial- and hnal-state mesons all carry 
0{1) fractions of the total soft and collinear momenta, respectively. Such conhgurations 
are represented by operators with fermion content Xc(... )XcQs(. • • )TCs. We will then consider 
“endpoint” conhgurations using the generalization in Table 2. These conhgurations occur when 
the momentum fraction carried by one of the partons tends to zero, so that the parton may be 
absorbed from the initial into the hnal state without hard momentum transfer. In particular, 
we will hnd conhgurations represented by operators with fermion content Xc(... )TCs- In both 
cases, using the counting rules for the building blocks containing soft-collinear helds described 
by the second column in Table 2, we will show that the operators representing the weak current 
at leading power do not contain soft-collinear modes. At leading power, soft-collinear modes 
appear only in time-ordered products of the weak current with subleading SCETn Lagrangian 
interactions, and with subleading terms in the interpolating currents for the meson states. 

The presence of the building block [n ■ dg n ■ dc)~^, which counts as an inverse power of A, is 
troubling at hrst sight. Naively, one could think that there would be inhnitely many operators 
of a given dimension and order in A. However, this is not the case: if an inverse derivative 
{n- dgfi- dc)~^ is added to a given operator, then it is necessary to also add two other building 
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blocks with d = 1 or one building block with d = 2 at the same time to obtain an operator 
of the same dimension. As is evident from the tables, this inevitably makes the resulting 
operator at least one power in A higher than the operator without the inverse derivative. In 
fact, from Table 1 we see that for operators involving only soft and collinear helds, with zero 
fermion number in both the soft and collinear sectors, the difference between the dimension of 
the SCETi operators and the order in SCETn power counting is given precisely by the number 
of occurrences of the building block [n ■ dgU ■ dc)~^- 

We focus hrst on the case of flavor non-singlet hnal states and will then discuss the modi- 
hcations necessary for the flavor-singlet case. We begin by considering SCETn operators with 
fermion held content Xc(... )Xc Qs(... )Tf 5 , corresponding to “typical” partonic conhgurations 
inside the initial- and hnal-state soft and collinear mesons. Using Table 1, and the fact that 
SCETi operators have dimension d > 3, it follows that leading-power contributions from these 
conhgurations are Later we will discuss other possible “endpoint” conhgurations, hnd- 

ing that they also appear at the same order in power counting. Erom Table 1 we see that 
with the exception of {n- dsfi- the building blocks satisfy [A] > d, so that leading-power 

operators of a given dimension must be generated with the minimal number of occurrences of 
this building block. Starting with the SCETi operator in (41) of dimension three, we hnd that 
the appropriate fermion held content cannot be obtained while remaining at d = 3 without at 
least one occurrence of the inverse derivative. The two possibilities at (!1(A^) are then 


r 


n ■ dr n ■ d. 


1 T ' 



n- dc 2 


n- ds 2 


{^fi. afi. <i} (46) 


and 




n ■ dr n ■ d. 


— o ^ 

n- dr 2 




(47) 


As in the tables, the notation is schematic: it is understood that the soft derivatives can act 
on any of the soft helds, and the collinear derivatives on any of the collinear helds. Using 
the equation of motion for the soft light-quark held, the above possibilities result in four 
independent operators, whose explicit forms are given in [12]. Their matrix elements can be 
expressed in terms of (endpoint divergent) convolution integrals involving twist-2 and twist-3, 
two- and three-particle LCDAs of the B meson and the light meson. The matching relations 
(46) and (47) are represented by the term shown at tree level in (26). 

Next, let us consider the SCETi current operators of dimension four. First, we observe 
that the operator qsAhciX"h does not match onto a leading order SCETn operator. Its soft 
bracket qsT"h is of order A^ and remains unchanged in the matching. The gluon held Ahc± 
must then match onto a d = 1 operator with collinear held content Xc(... )Xc. Inspection of 
the table shows that such an operator is of order A^, making the overall operator subleading. 
The only possibility to obtain a leading SCETn operator at d = 4 is 


to 

i 

1 T ' 

_ Ajp J. tAj n 


^ Qs^T'As 


n ■ dc 2 


_n-ds 2 


(48) 


At tree level, the matching (48) is represented by the term in (31). At dimension hve 

there are no possibilities for leading-power SCETn operators, due to the constraint [A] > d. 
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d 

[A] 


d 

[A] 

Xe 

3 

2 

1 

0 

3 

2 

2 

1 ^ 
n-dc 2 

-1 

0 

a 

3 

2 

5 

2 

Q. 

3 

2 

3 

2 


1 

3 

2 

1 f/- 

n-ds 2 

-1 

-1 

n-dcU-Dsc 

2 

2 

T" 

0 

0 

n-dsU-Dsc 

2 

2 


Table 2: Boost-invariant building blocks for SCETn operators containing non-zero fermion 
number and/or soft-collinear fields. The Dirac structures T" are defined in (43). a is the 
small-component projection of the soft-collinear fermion, as in (12). 


We thus need the SCETj operators only through d = 4 for leading-power matching. Finally, 
from the second column in Table 2, we note that replacing any of the soft or collinear helds 
in (46), (47), or (48) by soft-collinear fields results in power suppression. 

Our analysis has so far relied on the assumption that the held content of the SCETn 
operator is Xc(... )Xc □«(... )d{s, corresponding to “typical” parton conhgurations. We now 
consider possible “endpoint” contributions, corresponding to SCETn operators with fermion 
held content Xc(... )lKs. For this purpose, we consider the building blocks in the hrst column 
of Table 2, which generalize the hrst column of Table 1 to allow the possibility of non-zero 
fermion number in the soft and collinear sectors. Starting with the SCETj operator in (41) of 
dimension three, the leading SCETn operator is 

jA ~ _ ( 49 ) 


Again, from the second column in Table 2 we note that replacing any of the soft or collinear 
helds in (49) by soft-collinear helds results in power suppression. The operator in (49) can yield 
a leading-power contribution to the form-factor analogue of the correlator (5) when combined 
with leading-power meson currents and subleading Lagrangian interactions involving the soft- 
collinear modes [12]. Essentially, these interactions are summarized by the term 
in the ehective action of [25].^ Leading contributions can also arise from subleading meson 
currents and leading Lagrangian interactions. A contribution of this type to the B K*'y 
amplitude is illustrated in the hrst line of Figure 1. The interpolating current for the B meson 
takes the form 

TC,r"Q, + %sT"e + ‘ksV'a + ... 

(50) 




+ A^/2 + A' 


+ 


where the small-component projection a of the soft-collinear fermion held is related to 6 as in 

^More precisely, in the presence of the external weak current, the vacuum correlator of soft-collinear fields 
defining includes an extra soft-collinear Wilson loop Sl^Wsc [12, 32]. 
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(12). Similarly for the light meson, taking for example the pseudoscalar case, 


- i 




Xe|r'x, + 

+ 


+ 

A^ + 


(51) 


The subleading currents for both mesons suppress the contribution of XcT'dis by A^/^, so that 
it ends up being of the same order as the contribution of the four-quark operators. Finally, 
mixed cases can also occur, where an or 0{\), suppressed meson current from (50) or 

(51) is combined with an 0{X), or 0{X^^‘^), suppressed Lagrangian interaction, respectively. 
The relevant Lagrangian interactions in this case are given by O^+sc in [25]. By the 

same reasoning, we hud that all such endpoint conhgurations arising from the d = A SCETj 
operator in (41) are power suppressed. 

Before ending our discussion of heavy-to-light form factors, we consider the case of flavor- 
singlet hnal states. Operators corresponding to “typical” partonic conhgurations again have 
zero collinear fermion number, but may contain collinear gluon degrees of freedom in place of 
the fermion bilinear Xc(... )Xc. Requiring also that the collinear helds carry the appropriate 
twist and color quantum numbers to have overlap with the hnal-state collinear meson, there 
must be at least two such collinear gluon helds. From Table 1 we see that the new operators 
are obtained by the replacements 


1 W 
r 2 


{ct, 


AIc±-Ac_l{<9c^, 

n ■ dcTi ■ Ac Ac^ 


(52) 


in (46), and by the replacement 


1 T 

f’ J. c 


n-dc "2 




(53) 


in (47) and (48). From the leading SCETj current, there will also be new SCETn operators 
that combine with subleading soft-collinear Lagrangian interactions and meson currents to 
yield leading-power contributions. From Tables 1 and 2, we hud that at leading power the 
new operators are obtained by the replacement 

X, ^ Q^J—Iaci. (54) 

n- Os 2 

in (49). Although the right-hand side of (54) scales as X^A (compared to the left-hand side, 
which scales as A), leading contributions to form-factor matrix elements may still be obtained 
from subleading Lagrangian interactions involving the soft-collinear modes, which in this case 
are essentially summarized by the term of the effective action in [25]. 


3.2.2 General operators 

After this warm-up, we are ready to discuss the general case where the photon is not necessarily 
part of the SCETj operator. The new operators appearing in this case correspond to photon 
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emission from the spectator quark. The argumentation will be similar to the previous section; 
however, we will have to match also the n-hard-collinear part in (39): 

[“h-hard-collinear”] —x [“soft”] x [“soft-collinear”]. (55) 

The building blocks needed in this case are obtained from Tables 1 and 2 by exchanging n and 
h, dropping the collinear quark helds, and replacing the collinear gluon with the photon field. 
Note that the dehnition of the soft fields then involves Wilson lines in the h-direction. To 
distinguish them from the soft-helds appearing in conjunction with the n-collinear sector, we 
denote them by Tfj, Qj, and Ag. We also recall that the SCETi building blocks introduced in 
(33) are not invariant under soft gauge transformations; strictly gauge-invariant combinations 
are given by 

= -5l(a:-)X/,c(a:), Af/{x) = Sl{x_)AlXx)Ss{x _), (56) 

with the soft Wilson line dehned in (36). In general, the helds contained in the “n-hard- 
collinear” and “n-hard-collinear” brackets in (38) are not separately gauge-invariant. In order 
to match onto the building blocks in Tables 1 and 2 in the general case, we hrst translate to 
the gauge-invariant combinations appearing in (56). 

Let us again start by writing down a list of the relevant SCETi operators. By momentum 
conservation, they must have at least one n-hard-collinear and one n-hard-collinear held in 
addition to the heavy-quark held, and therefore start with dimension d = dn + dn (Td^) = 4: 


d = 4 = 


d = 5 = < 


3 + 1 : 

1 + 3 : 

4 + 1 : 

3 + 2 : 

2 + 3 : 
1+4: 
1 + 1 + 3 


XhcAjr^^rh (J^), 

(J^), 

(J^), 

Xhc ^ 5 


1 -p/ ‘P 

^.a—^hc 2 ^ ^hc 


XhcT'h, (57) 


j- _ ‘Y k p' w 


XhZ'^hcx di/icu r' h , 
ds Afici. •^~hc±_ r h. 


X^T'h (J^), Xj^^Ahc±Ahc±r h (J^), 


The symbols ..., in parentheses anticipate the notation to be introduced for these 
operators in Section 4. In constructing this list, we made the same simplihcations as in (41) 
for the form-factor case. In the above operators the held Aj^ stands for either the photon or 
a gluon held, which are treated on the same footing. In SCETn, we treat the photon as a 
collinear held in the n direction, ~ A. Since it appears only once in each operator, we 

are free to make such a scaling assignment. 

In the above list of operators, we have separately indicated the mass dimensions of helds 
in the n-hard-collinear, n-hard-collinear, and soft brackets, respectively. In those cases where 
the only soft held is the heavy-quark held, we have included it in one of the hard-collinear 
sectors in such a way that both hard-collinear brackets carry zero fermion number. SCETn 
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operators with zero fermion number can be constructed from the building blocks in Table 1 
which fulhll [A] > d. Beyond d = 5, operators appear which cannot be arranged to have zero 
fermion number in each hard-collinear sector. For example, 

d = 6; XhcT'Xj^^XhcT'h, Xj^J'Xhc Xj^J'h,... . (58) 


In Table 2 we have generalized the hrst column of Table 1 to include building blocks with non¬ 
zero fermion number, by splitting the various fermion bilinears in two halves in all possible 
boost-invariant ways. With the exception of Xc, all building blocks again satisfy [A] > d. In 
fact, since XcV'Xc = 0, the operator constructed from the n- and h-hard-collinear sectors must 
contain a factor (h ■ so that the bound [A] > d is recovered in the hnal operator. 

SCETj operators with d > 6 are therefore not relevant to a leading-power analysis. 

For the n-hard-collinear sector, we may use Table 1 to list the leading-order matching 
relations onto operators with minimal collinear held content Xc(... )Xc. This yields: 


dn — 1 • dl/jc_L 




d„ = 2 : 


1 y- ^■p/'Y' 


n ■ dcU ■ ds 


nr ' 


1 T 


n- dr 2 




n- dc 2 


A^ 


d„ = 3 : XhrT'h 


1 


n-dcTi-ds 


1 


n-dc n-da 


dr,. — 4 : X/icdl/,c±r'/i 


1 ‘Y 
n-dc 


1 T tv'r 


1 ‘X 

n-dc'^<^2^ 


ri'Cs 2 * 


Qsr"d{., 




n-dc 2 


n-da 2 


Q, ^F'FC, 


( 59 ) 


Note the presence of the soft Wilson lines, Ss {...) S']) ^r the dn = 1 case in (59). These 
factors are required in order to preserve soft gauge invariance, and can be derived via the held 
redehnitions (56). For the case of havor-singlet hnal states, we may again build additional 
operators using the replacements (52), (53). Also, in the cases dn = 1 and dn = 3, operators 
with collinear held content Aci_ appear at one order lower in A than those listed in (59), and 
can combine with subleading Lagrangian interactions to yield leading-power contributions, 
cf. (54). 

Similarly, in the h-hard-collinear sector, using the analogue of Table 1, we hnd the leading 
operators with h-collinear held content 


ds = 1 : 


dn = 2 ; 


•^hcL 


•^hcL-^hci. ) 


1 Y_kp'Y— 


A (em) 
•^c± 


A A : 

‘^c_L 


A 

A^ 
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dn 

dfL 


3 : 


, /l(em) 

i 

n ■ ds 

4 : 




-Q,^T"Ks 

a (e™) 




n ■ (9o 2 


dsi. 1 } 


A^ 

(60) 


Note that the soft Wilson lines appearing in the SCETn building blocks in (60) are in the 
opposite direction compared to those in (59). 

Returning now to (57), we hnd that dimension-four operators with d = dn + dn = ^ + ^ (J^) 
OT d = 1-1-3 {J^) can contribute at leading power. Similarly, at dimension Eve, those operators 
with d = 4 + 1 (J®) or d = 2 -h 3 {J^, J^) can contribute at leading power. The operators 
with d = 3-|-l, d = 4-|-l and d = 1 -|- 1 -|- 3 have been treated already in Section 3.2.1. They 
correspond to the case where the SCETi operator contains the photon held. The remaining 
operators, with d = 1-1-3 and d = 2-1-3, represent new contributions corresponding to emission 
of the photon from the spectator quark. 


4 Matching and factorization 

In the previous section, we have found all effective-theory operators that can contribute to the 
decay amplitude at leading power. Our analysis was concerned with the field content of the 
operators and the occurrence of inverse derivatives in SCETn, but we have not yet specified 
their color and Dirac structures. In this section, we present the relevant operators in all detail. 
We evaluate the Wilson coefficients necessary for the phenomenological discussion in Section 5 
and show that the resulting matrix elements can be brought into the form of the factorization 
theorem ( 1 ). 

4.1 SCETi matching 

We collect here the relevant SCETi operators as derived in Section 3. Again, we first consider 
the operators representing photon emission from one of the current quarks and then discuss 
those operators corresponding to emission from the spectator quark. We initially restrict 
our attention to flavor non-singlet final-state mesons. The additional operators that arise for 
flavor-singlet final states are considered separately at the end. 

4.1.1 Photon emission from the cnrrent qnarks 

Two SCETi operators are relevant for the case of photon emission from the current quarks, 
given by the d = 3 -|- 1 and d = 4 -|- 1 entries in (57). For the first of these, we write 

s, a) = 'Xhc{x+ + sn + Xi_) (1 -|- 75 ) {x_ + an + Xi_) /i( 0 ) . (61) 

In order not to overburden the notation, we refrain from indicating the flavor of the light- 
quark field. The dependence on the parameters s and a arises because the n-hard-collinear 
fields are allowed to live at arbitrary points on the h-light-cone, and the n-hard-collinear fields 
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at arbitrary points on the n-light-cone (cf. the discussion in Section 3.1). Furthermore, the 
position arguments of the helds have been multipole expanded, as appropriate for a product 
of fields (t)hc<hz<t^s- 


(t>hc{,x) (t)s{x) = (t)hc{,x+ + x±) (t>h^{x- + Xi_) 0^(0) + ... , (62) 


yielding the peculiar x dependence of the fields in (61). 

We use translational invariance to set a; = 0 and suppress the position argument in the 
following. The representation of the weak Hamiltonian for photon emission from the current 
quarks reads 


'wcurrent 


ds daC'^{s,a)J^{s,a) + 

j=l,2' 


ds / dr 


daCf{s,r,a)jf{s,r,a) + ..., (63) 


with the ellipsis denoting terms not relevant to a leading-power analysis. Here 

J^{s, a) = Xhc{sn){l + 75) .^^{an) h(0), 

Jf (s, r, a) = Xhc(sn)(l + 75) 4hc±(rn) h(0), 

Jf (s, r, a) = Xhc(sn)(l + 75) 4hc±(rn) ^(0) • 

We dehne Fourier-transformed Wilson coefficients as 

C-'{E,E^) = f<i‘> j iiae"" '’e“"''’’C"'(s,a), 

C^(E,E.„u) = jds jdr j Cf(s,r,a), 


(64) 


(65) 


where E = n-vn-P/2 and = n-P^/{2n-v). The quantity n-P is the large component of the 
total outgoing n-hard-collinear momentum, and similarly n ■ Hy is the large component of the 
outgoing photon momentum. We will suppress these quantities in the arguments of the Wilson 
coefficients in the following. The variable u denotes the fraction of the large component of the 
n-hard-collinear momentum carried by the quark held, and u = 1 — n is the fraction carried by 
the gluon held. The Wilson coefficients receive contributions from diherent weak-interaction 
operators, and we give separate matching results, and AjCf 2 , for the diherent Qi in 

(3). For b ^ s transitions we have 


^ps^pb C'i(hQCD) AfC'"^(pQCD, A'') + fr) ■ 

( 66 ) 

The same expression with s ^ d gives the coefficient for b ^ d transitions. Analogous expres¬ 
sions dehne AiC^ 2 - concentrate on the phenomenologically most relevant operators, 

which are Qj, Qi, and Qs- The scale /iqcD is the scale at which QCD and the ehective weak 
Hamiltonian are matched onto SCETi, and p is the renormalization scale in the ehective 
theory. 




Figure 4: Leading-order QCD diagrams for the matching of Qi and Qs onto . Other 
diagrams are power suppressed or vanish. 


The matching coefficients for Qi are obtained directly from the form-factor analysis and 
are given as 


ArC^ = 


erribE^ 


- 2 C ^, + -04 + 04 


ArOf 


erribE^ 

Stt^E 


erribE^ 

Stt^E 


Qcg+c®;) , 
(icg+cgj. 


(67) 


The tensor-current Wilson coefficients have been calculated through one-loop order, for C^j^(E) 
in [2, 33], and for C^^{E,u) in [33, 34], Explicit expressions for the combinations appearing 
in (67) are listed in Appendix A. In the above expressions, the MS quark mass must be 
evaluated at the QCD matching scale, i.e., nib = Eib{fJ,QCD)- For the process i? —> Dy, we have 
2E^ = 7715(1 — my/m%) and 2E = niB, with E^ and E dehned after (65). 

For the operators Qf (g = u, c) and Qg, we may deduce the one-loop matching onto A-type 
operators from results available in the literature [35, 36]. We hnd 



agCp 

All 


G^{x,)A^C^, 


AgC^ 


agCp 

Att 


GsAjC®, 


( 68 ) 


where Xg = frig/ml (we set = 0). The expressions for Gi{x) and Gg are the same as those 
in [18], and for convenience are reproduced in Appendix A. The 5-type matching is obtained 
from the diagrams in Figure 4, from which we hnd 

Figure 4(a): AJCf (u) = |l |, AJCf (u) =-A«Cf („), 

Figure 4(b): AgCf (u) = ^ ^ , A 8 C' 2 ^(u) = 0 . 

47 r^ 3 u 

The expression for f{x) is also given in Appendix A. 
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Figure 5: Leading-order QCD diagrams for the matching of Q 12 and Qs onto Jf. 


4.1.2 Photon emission from the spectator quark 

From Section 3, we also hnd leading-power SCETi operators corresponding to photon emission 
from the spectator quark. These contributions arise from dimension-two operators in the n- 
hard-collinear sector mapping onto purely collinear fields, cf. (59). The n-hard-collinear helds 
must therefore transform as a color singlet in order for the resulting operators to have non-zero 
matrix elements with the physical meson states. Also, for the chirality structure appearing in 
the Standard Model, only a single Dirac structure is relevant. Absorbing a factor 1/{2E) into 
the Wilson coefficients, and using the projection properties = h, = 0, the resulting 
four-quark operator takes the form 

Jf (s, r, a) = Xhc{sn){l + 75) ^Xhc{rn) X^(an)(1 -h 75)^^(0) • (70) 

In the presence of New Physics, additional operators can appear in the effective weak 
Hamiltonian. The class of non-standard operators includes four-quark operators with scalar, 
pseudoscalar, or tensor structures in place of the usual vector and axial-vector structures [37]. 
In addition to in (70), the following operators can then appear: 

7 = It=(l -75)^fc, 7 = 7=(1 +75)ffc. 

7 = 7o(l - 75)|Xio Xi;(l - 75)|fc , 7 = */.o(l + 75)7i|Xjo Xi;(l + , 

7 = X,.,(l - 75)7!;|xt„ Xjj(l - 75)|7x,.ft. (71) 

The representation of the weak Hamiltonian for spectator-quark photon emission is then 

Tf^ec'^ator ^ /ds f dr f da C^' {s, r, a) {s, T, a) + .. . . (72) 

k=i'^ 

In analogy with (65) for the H-type operators, it is convenient to introduce the Fourier- 
transformed coefficients 

C^{u) = jds jdr jda e^i^^+^r)n.P^ian-P, _ (73^ 
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The notation n ■ Pj = 2E^n ■ v anticipates that the h-hard-collinear quark held matches onto 
the photon (and a soft quark) in SCETn, see Figure 6. Clearly, Qj does not contribute to the 
matching onto C-type operators, and hence 

ArCf = 0 . (74) 


Evaluating the hrst two diagrams shown in Figure 5 for the operators Qf and Qs yields 


Figure 5(a): 
Figure 5(b): 


A\C^ 


N 471 [3 3 \4EE,^ ) 

Cp TTlh Us 1 
N 2E TT u 


(75) 


The function G(x, u) can be taken from [38] and is reproduced in Appendix A. For the charged 
decay mode B~ —> F“ 7 , the third diagram in Figure 5 also contributes: 


Figure 5(c): A’,‘Cp = 2 , 


(76) 


where q refers to the havor of the spectator quark inside the B meson. 


4.1.3 Flavor-singlet final states 

From (57) we hnd two new types of SCETj operators that can give rise to leading contribu¬ 
tions. For the chirality structure appearing in the Standard Model, the following operators 
are relevant: 


J^(s, a) = X^(an)(l -h 75 ) 4/^c±(s^^) h{0 ), 

J^(s, r, a) = X^(an)(l + 75) /i(0) (^f + ie^^) Al^p^{sn)Al^p^{rn ), 


(77) 


with and e'^ as dehned in (42). In writing the operator we have used the fact that 
at leading power the n-hard-collinear helds match onto purely collinear helds (and no soft 
helds) in SCETn, so that we may restrict attention to color-singlet operators in both the n- 
and h-hard-collinear sectors. Note that the relative sign of the e^^’^-term in the operator 
is without signihcance. The operator with the flipped sign is equivalent, if one also replaces 
C^{s,r,a) C^{r,s,a). Since the outgoing hadron is generated from gluonic degrees of 

freedom, the operators and can only contribute for flavor-singlet hnal-state hadrons. 

As usual, we dehne 


= I ds I da C^{s, a ), 


isn-P Aan-P^ / 


C^{u)= ds dr C^{s,r,a) 


(78) 


Qi does not contribute to the matching onto D- or F-type operators, and hence 

AjC^ = A-rC^ = 0 . 


(79) 


31 




The matching of the operator Qi onto and vanishes at zeroth order in and hence 


A\C^ = 0 , A\C^ = 0 . 


For the matching of Qs onto and we hnd 




Errib 


= 0 . 


(80) 


(81) 


4.2 SCETn matching 

We now write down the operators in the hnal effective theory and perform the matching of 
SCETi onto SCETn- The matching coefficients for this second step are called jet functions. 
We begin again with the flavor non-singlet case, considering photon emission from the current 
quarks as well as from the spectator quark. We then discuss the new ingredients needed for 
the treatment of decays with flavor-singlet hnal states. 


4.2.1 Photon emission from the current quarks 

The analysis in Section 3 showed which operator structures matches onto. The explicit 
form of these operators and their leading-order jet functions are given in [12]. However, since 
the non-factorizable part of the form factor can be simply dehned as the matrix element of 
the operator we do not need to perform this second matching step explicitly. 

The current operators Jj® and match onto 


Of{x = 0, s, t) = Xc{sn) (1 -h 75 ) ^Xc(O) Qs{tn) (1 - 75 ) ^IK^(O), 

(82) 

Of(a; = 0,s,t) = Xcisn) (1 + 75 ) | Xc(0) Qs{tn) (1 -h 75 ) , 

and two operators with color structure T“®T“, which have vanishing meson matrix elements. 
A consistent matching of SCETi onto SCETn beyond tree-level involves evanescent operators 
that mix with the operators in (82) [34]. Since we will be concerned primarily with an analysis 
at leading-order in RG-improved perturbation theory, and hence with matching coefficients 
only at tree-level, we do not list these operators here. The operators in (82) correspond to 
the d = 4 -|- 1 case in (59), (60). At tree level, the inverse derivatives appearing in (59) are 
accounted for via the relation 

1 r ^ 

_-_-_0(a;) =-? / ds (t){x + sn), (83) 

in - a + 


and similarly for n ^ n. Beyond tree level, the Wilson coefficients of the operators in (82) also 
develop logarithmic dependence on the light-cone variables s and t. As usual, we introduce 
the Fourier-transformed coefficient 


Di {uj,u) = / ds 




{s,t ). 


(84) 
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Figure 6 : Matching of Jf onto Of. Dashed lines denote hard-collinear fields, full lines soft 
helds. Note that the quark that emits the photon is described by a hard-collinear held in the 
h-direction, while the other two quarks are collinear with the n-direction. 


The Wilson coefficient of Of is a convolution of the SCETj Wilson coefficient Of with a jet 
function J7j_, 


Df{uj,u,fi) = 


00 


/■I / 2,Eoo 

'0 V F 


(85) 


The operator Of involves the jet function J\\. At tree level the two are identical, 


liTjl (ti, u)tree <j7d(^W)tree 


AtiCfOIs 


N 2Eu 


-5{u — v). 


( 86 ) 


The one-loop results for the two jet functions can be found in [34, 39]. We may recall that in the 
form-factor analysis the hard-scale matching coefficients are constant at tree level, independent 
of momentum fractions. In this case, up to hard-scale radiative corrections, expressions such 
as those appearing in (85) collapse into a simple integral over the jet function. Convolution 
with the meson LCDAs then yields a universal function Hm, identical for all form factors 
describing the same hnal-state meson M [39]. In contrast, for the B —I /7 analysis we see 
from (69) that even at tree level the coefficients are momentum-fraction dependent, so that 
the approximate universality represented by Hm cannot be utilized in this case. 


4.2.2 Photon emission from the spectator qnark 

For flavor non-singlet mesons, the only relevant SCETi operators not already present in the 
form-factor analysis are Jf. Eor the operator , the corresponding SCETn operator is 

Of(a; = 0,s,t) = (n ■ u)^ Xc(sh) (1-h 75) |Xc(0) Q5(th) (1-h 75) 4j!r^(0)|^5(0). (87) 

The remaining C-type operators corresponding to non-standard interactions are 

Of = (n ■ vf Xe (1 + 75) |xe Q.- (1 - 75) , 

03^ = (n ■ vf X, (1 - 75) |x, Q, (1 + 75) , (88) 

Of = {n ■ vf Xe (1 - 75) |x, Q, (1 - 75) , 
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= (n ■ vf Xe (1 + Qs- (1 + 75 ) , 

= (n ■ vf Xe (1 - 75)7^^|xe Q.- (1 - 75 ) • 

Note that the soft building blocks in (82) involve Wilson lines in the n-direction and a factor 
n/i next to Qs, while the soft helds in (87) and (88) involve Wilson lines in the h-direction and a 
factor ^ next to Q^. As a result, the matrix elements of the soft parts of the C-type operators 
will involve the same 5-meson distribution amplitude as the matrix element of . We dehne 
Fourier-transformed Wilson coefficients (recall that n ■ v = 1/n ■ v) 

Dficu, u) = jds jdt t ), (89) 

and 

{u, u, /i) = ^ (^n /i^ Cf {u, b ), (90) 

where Cq = 2e/3 for an up-type quark, and Cq = —e/3 for a down-type quark. From the 
Feynman rules of SCETi it follows that is proportional to the unit matrix, 

= (91) 

To see this, we recall that the n-hard-collinear and n-hard-collinear parts of the operators Jf 
match independently onto SCETn- The different operators Jf, and also Of, are distinguished 
only by the chirality of the fermion helds, and by the Dirac structure next to the heavy quark, 
both of which remain unchanged in the matching procedure. From the diagram in Figure 6, 
we then End 

J^' = l + OM . (92) 


4.2.3 Flavor-singlet final states 

Finally, let us discuss the new ingredients involved when havor-singlet final states are con¬ 
sidered. The modifications in this case are of two types. Firstly, new operators appear in 
the matching of A- and 5-type operators onto SCETn, corresponding to new contributions 
to form factors [40]. The new A-type operators are related to those appearing already in the 
havor non-singlet case by the replacements (52), (53) and (54) in (46), (47) and (49), respec¬ 
tively. Similarly, the new 5-type contributions are given by the replacement (53) in (48). 
The symmetry relations obeyed by the A-type form-factor contributions remain unchanged in 
the havor-singlet case; these contributions derive from SCETi currents = XhA'h, and the 
symmetry relations follow directly from the projection properties of the spinor helds Xhc and 
h. The new 5-type form-factor contributions are factorizable, involving the same 5-meson 
LCDA, and the leading-twist two-gluon LCDA of the light meson. We concentrate here on 
the second new ingredient in the havor-singlet case, namely the operators and . These 
operators contribute only to havor-singlet decays, and their contributions are unique to the 
radiative 5-decay mode. 
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Figure 7: Leading-order SCETj diagrams for the matching of the operators and onto 
SCETn operators and O^. Dashed lines denote hard-collinear helds, solid lines soft helds. 


The operator matches onto operators with collinear held content Xc(... )Xc: 


Of = (n ■ vf Xe(sh)(l + 75)|Xe(0) Q,(Fh)(l + 75)4(0)^s(0)4r^(0)4,^M|d{,(0), 

of = (n ■ vf tfsnfl - 75)|Xe(0) + 75)4(0)^.(0)47fo)4,^M|h{,(0), 

(93) 

and also onto operators with purely gluonic collinear held content, given by the replacement 
(53). The matching conditions in this case take the form 


^n/ ; X Co dvra^ , 2E^uj' \ ^n/ x 




(94) 


where we dehne 


Df{uj,uj',u) = ds dt 


(95) 


The jet functions, obtained from the Eeynman diagram in Figure 7(a), are given by 


The operators Of are of leading power despite the fact that their soft part involves an ad¬ 
ditional gluon held. The matrix elements of the corresponding operators involve non-valence 
Fock states of the B meson, but the presence of the extra gluon held is compensated by an 
additional inverse soft derivative. In a purely diagrammatic analysis, such a contribution can 
easily be missed, while the operator analysis performed in the previous section guarantees that 
all leading operators are included. As discussed after (59), also matches onto 

Of = {n-v) Qftnfl + 75)4(0)^s( 0)44(0)y±(0)y (0), (97) 
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which appears at one power in A lower than Of and Of. Similar to the A-type current in 
(49), this operator describes an “endpoint” conhguration of the two mesons where some of 
the partons carry very small momenta. It can combine with subleading terms in the SCETn 
Lagrangian or subleading meson current operators involving soft-collinear helds. The time- 
ordered product of Of with the subleading Lagrangians Cs+sc and Cc+sc contains terms with 
the same held content as Of and Of. The interactions with soft-collinear helds /Is+sc and 
Cc+sc [25] are both suppressed by A^/^ which makes the overall contribution leading power. Its 
presence signals an infrared divergence at a; 0 in (94), when the soft gluon in (93) becomes 
soft-collinear. The analysis based on power counting of soft-collinear modes provides a system¬ 
atic procedure to determine the presence or absence of such endpoint singularities. The iA-type 
contribution cannot be expressed in factorized form in terms of a hnite convolution integral 
over (generalized) meson LCDAs, and it is also not related to the non-perturbative quantities 
appearing in the form factor. Because these operators would contribute to B* P'j, we con¬ 
clude that for havor-singlet hnal states this decay mode does not obey a factorization formula 
such as (2). 

The operator is related to J*" by the replacement 

'Chc c-^r. C'Xhc —^ •A.hc±^hc± ■ (98) 

2n ■ dhc 

It maps onto a SCETn operator related to O^ in (87) by the corresponding replacement in 
(53). The jet function in this case arises from the Feynman diagram in Figure 7(b), and is 
identical to J'^ in (90). 

4.3 Matrix elements and factorization 

The analysis of Section 3 determined the SCETn operator structures onto which the SCETj 
current can be matched. Performing the matching explicitly and taking matrix elements 
yields expressions involving endpoint-divergent convolution integrals [12]. These infrared di¬ 
vergences indicate a sensitivity to endpoint momentum conhgurations, as verihed by the pres¬ 
ence of soft-collinear momentum regions at leading power. The contributions from such in¬ 
frared momentum modes spoil factorization and cannot be calculated perturbatively. Since 
we cannot reduce the A-type contribution to simpler hadronic quantities, we simply dehne the 
SCETi matrix elements® 

(M(p)|X;,erh|5W(n)) = -2EC,M{E)ii\MM{n)TMB+fv)\ , (99) 

where, as usual, 2E = n ■ v n ■ p. We have used the spinor wave-functions appropriate to the 
heavy-quark and large-energy limits: 

= = ( 100 ) 

®Since the quantity Cm(E )/is independent of nib, the heavy-quark flavor symmetry could be made 
manifest by extracting an additional factor in the definition of Cm{E) [41], similar to the definition of 

F{p) in (101) below. However, since we are concerned primarily with B mesons, we will use the normalization 
(99) that is commonly used in the literature. 
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where rj is the polarization vector in the case of vector mesons, and M. = 

The same power-counting arguments in Section 3 showed that subleading soft-collinear in¬ 
teractions are absent from the matrix elements of -B-type operators. Let us note, however, that 
the SCETn Lagrangian still contains leading-power interactions of soft-collinear gluons with 
both soft and collinear helds. In analogy to the decoupling of the soft-gluons in SCETi [31], 
it is possible to perform held redehnitions that remove the soft-collinear interactions from the 
leading-order soft and collinear Lagrangians [32], Under the same held redehnitions, the soft- 
collinear interactions also decouple from the operators Of and Of, and therefore the matrix 
elements of these operators factorize at leading power; the corresponding correlator diagrams 
consist of non-interacting soft and collinear parts (see Figure 1). The matrix elements of Of 
and of can then be written as convergent convolution integrals over the meson LCDAs: 


{0\Qs{tn) ^r7{s(0)|5(*)(t;)) = ^Tn^tr 




T J\4.b{*) (y) 




(M(p)|x,(sh)r|xe(o)|o) 


4 


n ■ ptr 


Muiji) V 


du /i). 


( 101 ) 

Let us note that the soft-collinear interactions are present, and do not decouple, for the 
operators with color structure T“ ® T“. The simple fact that an operator can be written as a 
product of soft and collinear helds does not guarantee factorization. 

We now collect the various elements and write down the leading-power decay amplitudes. 
For the Standard Model prediction we have 


X J dvj± V, In Uf(n,/i) 


= 2mB 

{V-fB\nw\B{v)) = Q, 
{P-1l\Hw\B\v)) = 2mB* 


Cvl + ^ f 4 >b ® ® ® Cf 


'lA^ I rrCnC 


c^Cp + 


-4>b ® fp(pp ® 


{Pir\'H„\B'(v)) =- ^,#.8 0 fpi.p ® Jj| ® C 2 ® 


( 102 ) 


The factors of 2 appearing in the above matrix elements arise from evaluating the polarization 
sums. For example, in the decay B U 7 the prefactor is = “2 if both the 

photon and the light meson have left-circular polarization, and zero otherwise. The metric and 
epsilon tensor in the transverse plane were introduced in (42). It is interesting to note that 
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for right-circular photon polarization, the B* P 7 amplitude is completely factorizable [42]. 
The same is not true for B* decays with left-handed photon polarization, and here spectator 
emission gives rise to isospin violation at leading power. 

For flavor-singlet hnal states, new non-factorizable contributions, which are not already 
present in the form factors, arise from the P-type operators. These operators contribute to 
the process B* P 7 for left-handed photon polarization, and we thus conclude that the 
amplitude for this process (in the flavor-singlet case) does not obey a factorization formula 
of the form (2). For right-handed photon polarization, on the other hand, the amplitude for 
B* —> P 7 remains factorizable in the flavor-singlet case. 

In the presence of New Physics operators with a chirality structure different from the 
Standard Model, one obtains additional contribntions 


{V^l\'H^\B{v)) = -mfFcj)B ® fvAv^ ® , 

{V^R\n^\B{v)) = +mfF^B ® fvAv^ ® . 


(103) 


A contribution to the non-standard Wilson coefficients C 5 and Cq would give a leading-power 
isospin-violating contribution to the P —> I /7 decay amplitude, with left- and right-handed 
photon polarization, respectively [15]. 

The factorization formulas for the B ^ V (and B* P) form factors at zero momentnm 
transfer involve the same hadronic parameters as appear in (102), bnt with different Wilson 
coefficients and . For example. 


{l/(p',7)| s7^6|P(p)) 


2i e^'^P^rilp'pP^ 


V{q^) 

rriB + my 


(104) 


2ie^‘'p°'p*p'^p 


niB 


C,7vi(B) + 


y/m^F 


<pB ® fx* 4>k* ® Ji. ® C 


where q = p — p'. Factorization theorems for all form factors are given in [39]. One can thus 
eliminate the non-factorizable piece ( 102 ) in favor of a form factor at = 0 , and rewrite 

the resulting expression in the form of the factorization theorem (1). Note that any choice of 
the renormalization scale p, will lead to large perturbative logarithms: the coefficients and 
contain logarithms of the hard scale, ln(/i^/m^), and the jet fnnction contains logarithms 
of the hard-collinear scale, /umh)■ One can resnm these logarithms by solving the RG 

equations for the Wilson coefficients and the jet fnnctions [39]. The phenomenological impact 
of this resnmmation will be discussed in Section 5. 


4.4 Light-quark masses 

We have demonstrated factorization for P —R 7 by expanding the weak Hamiltonian onto a 
complete basis of operators in the effective theory, and then isolating those contributions that 
cannot be absorbed into the B ^ V form factor. For these contribntions, we demonstrated 
the insensitivity to infrared momentum regions that would signal endpoint divergences in the 
hard-scattering convolution integrals. This infrared insensitivity in tnrn was demonstrated 
by the deconpling of the soft-collinear messenger modes that could potentially communicate 
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between the soft and collinear sectors of the theory to spoil factorization. An interesting 
question to ask is how light-quark mass terms can affect the conclusions drawn from analyzing 
the massless case. 

Light quark masses of order A can be expanded in the propagators of hard-collinear par¬ 
ticles, 


so that the quark masses can be ignored at leading power for such momentum regions. In 
contrast, the propagators of collinear or soft particles. 


pl-ml 


pl-ml 


(106) 


cannot be expanded, so that the quark masses appear at leading power in the low-energy 
theory. For the region of soft-collinear momentum we have 


1 


P 


2 

SC 


— rrig 



(l + , 


(107) 


so that in the presence of such a light-quark mass, no pinch singularities arise in the fermion 
propagators from this region. The absence of the soft-collinear mode does not imply that all 
quantities in the low-energy theory factorize, but rather that the question of factorization has 
become more subtle. Individual Feynman diagrams for the soft and collinear modes contain 
divergences that are no longer regulated by dimensional regularization. Additional regulators 
may be introduced to make the diagrams individually well-dehned, but such regulators link 
the soft and collinear sectors. Demonstrating factorization then involves showing insensitivity 
to the additional regulator. 

We restrict ourselves here to the more tractable case where rUg A. We recall that the 
fundamental object under study, the correlator (5), is free of infrared singularities, and so has 
a smooth limit as rUg 0, keeping A hxed. The dependence is not analytic, however, so that 
we cannot simply treat the mass term as a perturbation, and work systematically to arbitrary 
order. The non-analyticity is associated with new regions that appear in the presence of 
the mass term. Such non-analyticities arise only from the propagator denominators. In one- 
loop examples, we hnd that the leading nonanalyticities are quadratic in the fermion masses 
as long as rUg A. Terms linear in the mass appear only from the numerator structure, 
with denominators described by the massless case. Assuming that this property persists 
at higher order in perturbation theory, it is straightforward to show that the soft-collinear 
helds decouple in the usual way from the new terms in the soft and collinear Lagrangians. 
Factorization properties of the decay amplitudes are therefore unchanged, except that the 
hadronic parameters (Cv^, 0 b) are modihed by quark-mass effects. Beyond linear order, 

the analysis may become more complicated. If the leading non-analyticities are quadratic in 
rUg, the factorization formula (1) holds at least up to terms of 0{rr?J 
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5 Phenomenology oi B —^ 1^*7 


From the decay amplitude, we obtain the following result for the B K*'j branching fraction: 

Tsms 


Bt{B K*'^) = 


An 




m 


B 


i^r, 


where we introduce the notation 


A = ^ ^ ® fKi4>Ki ® >J± ® C'f 

= -^soft “1“ -^hard • 

Neglecting contributions proportional to Vub, the Wilson coefficients are 

CpOis ( 


Gf. 




^7 + 


An V 


GsGs + CiGi(xc) 


cf = 


Gf. 




cb 




ml 


1 / J 


^^C■' + 0(al) , 
ArCf + 0(a,). 


( 108 ) 


(109) 


( 110 ) 

( 111 ) 


In the above equations, we have suppressed the scale dependence of the various quantities. 
The Wilson coefficients of the effective weak Hamiltonian, Ci, Cy, and Cg, depend on the 
renormalization scale /iqcD in the “full theory” consisting of ordinary QCD and the operators 
Qi. The quantity depends on /Xqcd as well as on the renormalization scale in SCETi, 

/i. This dependence on /i is canceled by the opposite dependence of Cic*(At)- Since we will 
determine the non-factorizable part .Asoft directly from a physical form factor, it is simplest 
to choose the scale /iqcD = ~ rufe in this part. This choice guarantees the absence of large 

logarithms in the perturbative expansion of G"^. Such a choice is not appropriate for the 
factorizable part, Alhard, but since the two parts are separately RG invariant, we are free to 
choose different scales in the two parts [39]. 

Since we have performed two matching steps for the factorizable part, a single choice 
for all renormalization scales inevitably leads to large perturbative logarithms: the SCETi 
matching coefficients depend on the hard scale mb, and the jet function on the hard-collinear 
scale y/ojmb ~ y/Amb ~ 1.5 GeV. Also, the meson LGDAs are typically given at a low 
renormalization scale /i ~ 1 GeV. By solving the RG equation for the Wilson coefficients, 
one can sum up perturbative logarithms of ratios of these scales, and match consistently onto 
the hadronic matrix elements at the low scale. Below we will give the result obtained after 
resummation and compare it to the hxed-order result. 


5.1 Non-perturbative input 

In order to evaluate the branching ratio, we need the value of Cic* at the kinematical point 
= 0, corresponding to maximum recoil energy of the K* meson, as well as the meson LGDAs 
05 and (j)K*i to evaluate the factorizable part. Unfortunately, there is no direct experimental 
information on these quantities available, so that we will rely on sum-rule determinations. The 
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value of Ce'* can be determined from any B K*^ form factor at zero momentum transfer, 
because all such form factors have the same non-factorizable part. We will use the vector 
form factor V, which fulhlls the factorization theorem (104). This choice is convenient, since 
the factorizable part of V is 0{as) suppressed compared to other form factors, e.g. the tensor 
form factor Ti. In other words, Cy vanishes at tree level, and we have 

_m^ ^ CKi{E,fi) + 0{as{mh)as{^/]^b)) • ( 112 ) 

tub + my 

The B K* form factors have been determined from light-cone sum rules. The most recent 
evaluation gives 1^(0) = 0.411 ±0.045 [43]. This value is compatible with, but somewhat lower 
than the earlier result V^(0) = 0.458 ±0.069 [44]. Using the known one-loop expression for Cy 
[2, 33], the sum-rule determination yields 

a* =Cxi(^= ^,h = m„) =0.40±0.04. (113) 

Let us check whether the sum-rule results for the axial form factor Ai and the tensor form 
factor Ti give the same value of C,k*^, as required for consistency with the heavy-quark limit. 
The relations between the vector form factor V and the axial form factor Ai is especially 
simple, i.e.. 


ml U(0) 

{ms ± mvY ^i(O) 


1 ± 0{l/mb ). 


(114) 


This form-factor relation does not receive perturbative corrections [39]. The sum-rule value 
for this ratio is very close to unity, and the value Cx* = 0.39 ± 0.05 extracted from the axial 
form factor is consistent with the value obtained from V. The relation between the vector and 
tensor form factors is slightly more complicated, since their Wilson coefficients and 
are different. In particular, the factorizable piece is not suppressed in the case of the tensor 
form factor, for which the sum rule evaluation gives Ti(0) = 0.33 ± 0.04 [43]. Evaluating the 
factorizable part using the hadronic input as given in Table 3, and including the resummation 
effects as discussed in Section 5.2, we obtain = 0.37 ± 0.04, again consistent with the 
value from U(0). 

The iL*-meson LCDA is used as an input for the sum-rule evaluation of U(0). Following 
[43], we parameterize this function in terms of the lowest two Gegenbauer moments as 


(f)Ki{u,fi) = 6 m (1 


u) 


1 ±ai(/i)C'i/^(2u 


l)±a2(/i)C'2/"(2u-l) , 


(115) 


and we use ai(l GeV) = 0.1 ± 0.1 [45], 02(1 GeV) = 0.1 ± 0.1. For the K* decay constant we 
use /ii'*(lGeV) = 170 ± lOMeV. For the B meson, we take the model [46] (see [47] for an 
alternative form) 


Ai = 1 GeV) 




UfJ, 


TT + fJ? \_UJ^ ± /i 




2((Tb - 1) UJ 

-s 


(116) 
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rriB 

5.28 GeV 

Tb = {tb+ ± tbo)/2 

1.60 ps 

mb{mb) 

4.25 ±0.1 GeV 

mb 

4.8 GeV 

Ib 

200 ± 30 MeV 

Cki 

0.41 ±0.04 

^B 

460 ± 110 MeV 

(Tb 

1.4 ±0.4 

mx* 

894 MeV 

fK* (1 GeV) 

170 ± 10 MeV 

ai(lGeV) 

0.1 ±0.1 

02(1 GeV) 

0.1 ±0.1 

^'■QCD 

217 MeV 


0.040 ±0.002 

rric 

1.1 ±0.2 GeV 

Climb) 

1.108 (LL) 


-0.320 (LL), -0.311 (NLL) 

Cf^\mb) 

-0.151 (LL) 


Table 3: Numerical input values and uncertainties. See [50] for the definition of Leading- 
log (LL) accuracy is sufficient for Ci and while we need the next-to-leading-log (NLL) 

value for in C^. 

with parameters \b = 460 ± 110 MeV and as = 1-4 ± 0.4. To leading order in perturbation 
theory, and if no RG improvement is performed, the factorizable part depends only on the 
hrst inverse moment of the LCDA, dehned as 

—0B(o;,/i). (117) 

Jo ^ 

Note that in the model (116), \b{i^ = 1 GeV) = \b- The quantity F dehned in (101) is related 
to the R-meson decay constant, ydrzB/s = KF{fJ')F{fj,), up to higher orders in l/rrib, with [48] 

KF{^x) = l + ^^^J^(3\n—-2]. (118) 

Att \ J 

We use the value = 200 ± 30 MeV for the R-meson decay constant, which lies in the 
ball park of lattice and sum-rule determinations of this quantity. For the &-quark mass in 
the MS scheme, we use mb{mb) = 4.25 ± 0.1 GeV [49]. For the charm-quark mass we use 
rUc = 1.1 ± 0.2 GeV, a range of values that corresponds to an MS-mass with a scale between 
rric and mb- 

Note that the branching ratio depends only weakly on the value of the pole mass mb = 
4.8 GeV, through the one-loop corrections to C^. The pole mass can be eliminated in favor 
of a low-scale subtracted 5-quark mass. For the R-meson lifetime, we use tb = 1.60 ps. We 
use three-loop running for as with Aq^]], = 217 MeV [49] and work with n/ = 4 below the 
scale 11 = mb and n/ = 3 below the intermediate scale fi = \fKbm^ = 1.55 GeV ~ rric, where 
A/i = 0.5 GeV represents a typical hadronic scale. Numerically, using n/ = 3orn/ = 4 makes 
very little difference in the results. 

5.2 Resummation 

The anomalous dimensions of the SGETj cnrrent operators were calculated in [39]. The 
anomalous dimensions of the SGETn fonr-quark operators are given by the anomalous dimen- 
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sions of the meson LCDAs. In both SCETj and SCETn, the operators are non-local along one 
or more light-cone directions; their anomalous dimensions are distributions that describe how 
the operators at different light-cone coordinates mix among themselves. The necessary steps 
to solve the evolution equations were spelled out in detail in [39] and we refrain from repeating 
the discussion here. After resummation, the factorizable part of the amplitude takes the form 

Aard = 0^(/^)<g> A* (h)0xy (jU) Hh) (/ife), (119) 

where Ui evolves the coefficient from the hard scale to the intermediate, hard-collinear 
scale, and Un describes the second evolution step down to the low scale of order 1 GeV. Aard 
is independent of the scales /Xqcd, At/i, and /x, up to higher orders in as at these scales when 
evaluated at hxed order in perturbation theory. To estimate the uncertainty from higher-order 
perturbative contributions, we will independently vary the hard and the hard-collinear scales 
by a factor of \/2 around their central values and fXi = y/AhUib with Ah = 0.5 GeV. 

Throughout, we set /iqcD = /U/i- Let us stress again that we count nic as a hard scale and 
therefore do not resum perturbative logarithms of rric/mb- In view of the fact that rric is 
numerically rather close to the intermediate scale it might be advantageous to count rric as 
hard-collinear instead of hard, which would allow one to also resum such logarithms. However, 
an effective-theory framework for such a treatment is not yet available. 

To gauge the size of the corrections to the amplitude, we will express them in terms of the 
leading-order result 

^(0) ^ (jA{0) ^ ^ E^mbinib) C^^{mb) A* • (120) 

Including 0{as) corrections, and normalizing with respect to the leading-order result, the 
contribution of the non-factorizable part is 

^ ^ ^(0) ^ Q_Q52 ^ 0.027) + i(0.062 ^ 0.014 ^ 0.016)] , (121) 

where the hrst uncertainty comes from varying /Xqcd = Ah and the second from the variation 
of rric in the contribution from Qi. The larger values of Aoft arise from the lower values of 
the scale /i/i, which we indicate with the symbol “=f” . The factorizable part of the amplitude 
is smaller, similar in size to the 0{as) correction to the non-factorizable part. We hnd 

= A°^ [(0.055 T 0.010 T 0.009 ± 0.019) + f(0.031 ^ 0.005 ^ 0.004 ± 0.011)]. (122) 

The hrst uncertainty comes from varying the intermediate scale /i*, the second from varying the 
hard scale fih- The third uncertainty is estimated by varying the hadronic input parameters 
within the ranges in Table 3. It is dominated by the uncertainty in the H-meson LGDA and 
in rric/nib, see Table 4. 

Let us compare this result to what is obtained in a hxed-order calculation. Using a common 
scale Ah = Ai = V^hnib, we hnd 

Aard = ^^°U0.116 + ffi.062). (123) 
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fs 

As 

CTb 

Iki 

ai 

0,2 

rriclmb 

103AAard/Al(°) 

±8±5i 

=F14 =f Si 

TSTi 

±3±2i 

±5±i 

±5 =F i 

5-7i 

-10+5i 


Table 4: Hadronic uncertainties in the evaluation of the factorizable part .Ahard in units of 
10“^ A.The parameters are varied within the bounds given in Table 3. We add the sym¬ 
metrized errors in quadrature to estimate the total hadronic uncertainty for which we obtain 
AAard = ^(°H0.019 ± O.OIH). 


The above scale choice guarantees the absence of large logarithms in the jet function, but 
will lead to perturbative logarithms ln(/i^/m^) in C^. We could instead use a large scale 
fJ'h = f^i = uif,, which eliminates the logarithms from but induces logarithms AhTUb) 

in the jet function. This gives 


Aard = .4<“> (0.032 + iO.019). 


(124) 


With this scale choice, the result is more than a factor three smaller. The hxed-order cal¬ 
culation suffers from large scale uncertainties, which are greatly reduced by performing RG 
improvement. 

The one-loop corrections to the jet function are beyond the accuracy of our calculation, 
since at the same order also the effect of the two-loop running and the 0{as) corrections to 
would need to be included. However, to get an idea of the impact of these corrections, we 
evaluate Aard including the one-loop jet function on top of the leading-order evolution. The 
effect of perturbative corrections to the jet function are moderate: 

= dt*"’ 1(0.063 T 0.007 T 0.010 ± 0.023) + 1(0.034 q: 0.003 T 0.004 ± 0.013)], 

(125) 

with the errors arising from the same sources as in (122). The uncertainty from varying 
the renormalization scales fih and fit in the factorizable part is very small, indicating that 
perturbative corrections to the above result are likely to be small. 

Combining (121) and (125), we obtain from (108) the branching ratio 


Br(R ^ Ay) 


6.6 ± 1.3^^^ ± 1.3A/m6 0.7 ckm ± 0.7,,.,^^ ± 0.4 


X 10“®. 


(126) 


We have isolated the uncertainties associated with Ay i power corrections, the CKM prefactor, 
and the renormalization scale; the latter uncertainty is dominated by the soft term (121). 
The final uncertainty in (126) is associated with the remaining input parameters. If is 
determined from the tensor form factor, as discussed after (113), the branching ratio comes 
out lower, with central value 5.7 x 10“^. 

The uncertainties from power corrections are difficult to quantify. We have followed stan¬ 
dard practice and estimated these corrections by assuming A ^ 500 MeV, which gives a 10% 
uncertainty in the amplitude. In the case of charmless two-body decays classes of enhanced 
power corrections have been identihed which are larger than this naive estimate [16, 38]. How¬ 
ever, since these chirally enhanced and annihilation contributions are absent in our case we 
believe the power corrections are of natural size. In addition to the scale mb, our results 
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involve the hard-collinear scale iihc ~ V and the charm quark mass. The expansion of the 
amplitude in the inverse of the hard-collinear scale is quadratic, so that these contributions 
are of the same size as the A/rrif, corrections. The situation is similar for the power corrections 
associated with the charm quark. We have analyzed the diagram shown in Figure 4a and hnd 
that the corrections are Similar conclusions were reached in [18] based on the endpoint 

behavior of light-cone wavefunctions. While we do not have a formal proof that linear terms 
are always absent, this seems plausible since we are integrating out a heavy charm quark. At 
leading power in A/mb, the charm quark corrections are furthermore suppressed by ag at the 
hard or hard-collinear scale, which compensates for the fact that the Wilson coefficient Ci 
multiplying these contributions is approximately 3 times larger than the dipole coefficient Cy. 

Although we hnd agreement with previous analytic results in the literature [17, 18], there 
are some differences in the evaluation of the dual branching fraction (126) due to hadronic 
input parameters and to our RG analysis. The largest difference comes from the overall 
normalization given by the input value of the form factor. The remaining difference in the 
soft contributions from Qi and Qg is due to our smaller value of the charm-quark mass, and 
to our use of for the soft matrix element in place of Ti, which includes (higher-order) 
hard-scattering terms unrelated to Qi and Qg. For the hard contributions, the difference is 
accounted for by ~ la variations in the input values of nic, \b and Oi, and by our inclusion 
of a complete leading-order RG analysis. ® 

The B 77*7 branching ratio has been accurately measured by the GLEO [51], Belle 
[52], and BaBar [53] collaborations. An average of their results gives Br(R° —= 
(4.03 ± 0.26) X 10“^ and Br(R+ —i7*+7) = (4.01 ± 0.20) x 10“^. At leading power, the 
factorization formula predicts that the branching ratios for the charged and neutral decay are 
identical. At subleading power, photon emission from the spectator quark breaks the isospin 
symmetry. This effect was estimated in [15]. 

Our result for the decay rate is 65% larger than the experimental result, or 1.2 a with 
the errors in (126). If we were to take such a discrepancy seriously, it would be difficult 
to attribute the difference in the exclusive decay to New Physics, given that the prediction 
for the inclusive b ^ s'j decay agrees well with the experimental result [54]. In principle, 
it is possible that New Physics affects the two decay modes differently: spectator emission 
is suppressed by 1/ml inclusive decay, while it can be leading order in the exclusive 

decay if New Physics is present. However, the presence of such operators would typically 
lead to large isospin asymmetries, as we discuss below in Section 5.3. It is also not plausible 
that higher-order perturbative effects could account for the difference. Either the sum-rule 
result for Cx* is ~ 30% too large, or there are power corrections of this size which violate 
the factorization theorem (or some combination of these possibilities). If we instead use the 
experimental result for the branching fraction to determine the non-factorizable part of the 

®For example, in place of the quantity a^{K*j) in equation (55) of [18], we have 

a^{K*-f) -0.320 -h 0.009 -f [C/c* /Ti(0)](-0.098 - 0.023z) -h (0.055 + 0.031z - 0.042). 

The first two terms arise from the dipole coefficients and The third term represents the soft 

contribution of Qi and Qsj from (121), and similarly the fourth term gives the hard contribution of Qi and 
Qs, from (122). 
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form factor, we obtain C,k*^ = 0.31 ± 0.02. 


5.3 Isospin and CP asymmetries from New Physics 


We now consider the effect of the C-type operators, which arise from spectator emission. As 
discussed earlier, the two operators and Oq that contribute to i? —> K*j have vanish¬ 
ing Wilson coefficients in the Standard Model. Consequently, spectator emission is power 
suppressed in the Standard Model, and the New Physics effects associated with these oper¬ 
ators can lead to isospin and CP asymmetries that are enhanced over the Standard Model 
predictions. 

Because the jet-functions of the C-type operators are proportional to the charge of the 
spectator quark, the presence of these operators induces an asymmetry between the charged 
and neutral decay modes: 


P(5+ ^ K*-^) - r(B'^ K*°j) 




B 


IGttAI 




yd\2\ 


1.5 X 10 
tb 


-3 


l4“l" + bet - \\Cl? - jiqi" 1 X TeV‘ 


(127) 


In general, the Wilson coefficients are functions of the light-cone momentum fraction 
of the light current quark and are convoluted with the LCDA of the light meson: C = 
Jq du(l)K*^{u)C{u). If the New Physics takes the form of four-quark operators at the hard 
scale, the coefficients will be constant and If the New Physics is isospin symmetric, 

then = Cf and the difference (127) is positive. However, in general the operators entering 
the charged and neutral decays can have different Wilson coefficients. Experimentally, the 
difference between the decay rates, normalized by the average lifetime tb as in (127), is 
(—3 ± 3) X 10“®, where we take a weighted average of the values from [51], [52] and [53]. The 
most recent estimate of the difference in the Standard Model is (—5 ± 3) x 10“® [23]. From 
(127), we hnd that the Wilson coefficients evaluated at renormalization scale fi ~ mi, must 
obey |C* 5 , 6 | ^ 5 x 10“^TeV~^, meaning that even with present precision one is able to probe 
New Physics effects at scales of several TeV. 

Let us now turn to the time-dependent CP asymmetry. The and B^ mesons can decay 
into {K'l 'Jb) or These two states can be combined into a CP-even and a CP-odd hnal 

state, and since the spins are not measured, the experiments give the sum of the two rates. 
Expressed in terms of the amplitudes Al,r = A{B^ zj) and Al,r = A{B^ zj), 

the time-dependent CP asymmetry is 


T{B^{t) K*-f) - T{B^{t) 77*7) 

P(H0(t) ^ 77*7) + r(H0(t) ^ 77*7) 


( 128 ) 


[\Al\ 

\^+\Ar\^-\Al\ 

P — Mijp] cos(Am_Bt) — 2Im 

^AlAI + ArA*j,) 

[p J 

sin(AmBt) 

\Al\ 

\^ + \Ar\^ + \Al\ 

|2 + 
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The coefficients p and q relate the mass to the flavor eigenstates: \Bhl) = p\B^) ± In 

deriving the above expression, we have assumed \q/p\ = 1. This holds to good approximation, 
since the width difference in the B^ system is very small. 

In the Standard Model Al and Ar vanish to leading power in l/rrih. The coefficient of 
sin(AmBt) is thus power suppressed. The prefactor of cos(AmBt) also happens to be small in 
the Standard Model. The reason is that, up to terms which are doubly Cabibbo suppressed, 
the 6 —> s amplitude has only a single weak phase, so that there is no CP violation in the 
decay. The direct CP asymmetry from the Cabibbo suppressed terms is 

- UlP ^ (V^\ ( (/f-71 CiQ? + C^Ql + SL t , CXA* 1 

\ArP+\AlP \{K’~,\C,Q'i + C2Qi + El3QQi\B«)) ° 

(129) 

= ri\l [-0.14 ± 0.03a/„ ± 0.03„,.a. ± 0.04„, ± 0.02] + ©(Aj,). 

Power corrections are estimated as 20% of the leading result. Uncertainties associated with 
scale variation and the charm-quark mass are indicated explicitly, and the final uncertainty 
is due to the remaining input parameters. At leading power, the asymmetry is identical for 
neutral and charged B mesons. Here Ac = Vus ~ 0.22 is the sine of the Cabibbo angle, 
and the Wolfenstein parameter p is related to the imaginary part of Vub- The CP violation 
in the Standard Model is thus negligible. New Physics can change these predictions rather 
dramatically. If the New Physics has operators that induce leading-order contributions to the 
decay amplitudes Al and Ar, large CP-violation effects from the interference of mixing and 
decay can be observed, even if the New Physics operators do not have a new CP-violating 
phase [55]. In the presence of New Physics operators with additional phases, asymmetries in 
the decay can also occur. Note that the operators and Oq can contribute in both cases. 
These operators are suppressed by l/mf in the inclusive decay. A difference in the exclusive 
and inclusive direct CP asymmetries could be explained by the presence of such operators. 

Recently, both BaBar [56] and Belle [57] have performed measurements of the time- 
dependent CP asymmetry in the B K*'y decay. Within large errors, their results are 
consistent with a vanishing CP asymmetry. 

6 Discussion and conclusions 

We have established the factorization formula (1) for R —> Uy decays, which provides the 
basis for the phenomenological analysis of the decays B K*'p and B py. To perform the 
diagrammatic analysis of the factorization properties of the amplitude, we studied a current 
correlation function from which the B —> K*'p amplitude can be extracted. The diagrams 
contributing to the correlator are expanded around the heavy-quark limit using the strategy 
of regions. The different momentum regions are represented by corresponding fields in the 
effective theory. Similar to the case of two-body decays, such as B ^ nn, the decay B K*p 
involves energetic partons propagating in two directions. It is then necessary to introduce 
collinear fields along the light-meson direction as well as the photon direction, making the 
effective theory analysis more involved than in the case of the heavy-to-light form factors. 
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A large number of possible operator structures appear and it becomes crucial to have an 
efficient way of identifying the relevant operators both in SCETi and SCETn. There are three 
complications that make the construction of the operators non-trivial: (i) the power counting 
in the two effective theories is different and the power of an operator in SCETn does not impose 
a strong constraint on the SCETi operators that match onto it; (ii) the SCETn operators can 
contain inverse soft derivatives which count as inverse powers of the expansion parameter, and 
(iii) two collinear sectors are present. We have set up an efficient formalism to construct the 
operators which addresses all three issues. We classify the SCETi operators by their dimension 
instead of their power in A and derive a constraint on the maximum dimension of the SCETi 
operators beyond which they cannot match onto leading SCETn operators. We construct the 
SCETn operators from gauge-covariant and boost-invariant building blocks. The use of these 
building blocks makes it simple to identify how many inverse soft derivatives can occur. Finally, 
we separately match the collinear helds from the two sectors to account for the structure of the 
SCETi Lagrangian. Once the operators are identihed, it becomes possible to make all-orders 
statements concerning factorization by identifying which classes of operators are insensitive to 
infrared momentum regions. Our analysis shows that the B — K*'y amplitude indeed takes 
the form of the factorization formula ( 1 ). 

The basis of effective-theory operators which we constructed for the factorization proof 
can be used to analyze the effects of New Physics in the decay i? —I/ 7 . These operators 
also describe decays with a flavor-singlet hnal-state meson, and the related decay processes 
B* —> P 7 of B* vector mesons. The New Physics operators yield calculable contributions at 
leading power to the isospin asymmetry and the time-dependent CP asymmetry in P —^ -^* 7 , 
so that measurements of these asymmetries provide useful constraints on the corresponding 
New Physics operators. For flavor non-singlet hnal-state mesons, we found that the B* — 
P 7 amplitudes obey the generalized factorization formula (2). For the havor-singlet case, 
however, the P* —P 7 decay amplitude does not obey a factorization formula of the type 
(2), highlighting the non-trivial nature of the P —E 7 factorization formula. More precisely, 
only the P* — P'yi mode violates the factorization formula. The power counting of soft- 
collinear helds makes apparent the sensitivity to infrared momentum regions at leading power, 
and consequently the breakdown of factorization in this case. In contrast, the B* P'Jr 
mode is completely factorizable, involving only the second term of ( 2 ), and the same power¬ 
counting arguments show that the entire amplitude in this case is perturbatively calculable as 
a convergent convolution integral over meson LCDAs. 

The ehective theory approach allows us to disentangle the diherent scales in these decay 
processes and to resum perturbative logarithms of these scales by solving the RG equations 
in the ehective theory. We reanalyzed the P — K*'y branching ratio with recent values for 
hadronic input parameters taken from light-cone sum rules, and presented the hrst analysis 
of the hard-scattering terms to leading-order in RG-improved perturbation theory. We also 
evaluated the impact of perturbative corrections to the jet function, which are potentially 
signihcant because they arise at a low scale p ~ 1.5 GeV. Parts of these corrections were 
identihed in the diagrammatic analysis of [19], which isolated large logarithms occurring in 
the hard-scattering kernels. The result for these logarithms is precisely reproduced by inte¬ 
grating the one-loop jet function J7j_ in [34, 39] over the tree-level hard-scale coefficient. In 
our analysis the leading logarithms are resummed and the remaining one-loop jet-function 
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correction increases the factorizable part of the amplitude by approximately 15%. Depending 
on which QCD form factor is used to determine the SCET quantity Ck* , our result for the 
B — K*'^ branching fraction differs from the experimental value by la — 2a. The discrepancy 
is smaller than in earlier evaluations, mostly because the new sum-rule values for the B K*^ 
form factors are somewhat smaller than earlier results. 

The jet function appearing in i? ^ Dy is identical to that appearing in i? —> form 
factors. The universality of the jet function gives rise to new symmetry relations between 
different form factors when perturbative corrections at the hard scale are neglected. The 
reason is that the SCETj Wilson coefficients for the form factors are constant at tree level. 
In this approximation, the integral of the jet function over meson LCDAs yields a universal 
quantity, identical for all form factors describing the same hnal-state meson [39]. The same 
integral over the jet-function also appears in processes such as B ^ tttt at the same level 
of approximation [58]. In contrast, for B ^ V'j the tree-level hard-scale coefficient is not a 
constant. Here, even in this leading-order approximation the universal hadronic parameters 
must be taken as the meson LCDAs, and the perturbative expansion of the jet function is 
essential to retain predictive power in the factorization formula. 

In summary, we have presented the hrst factorization analysis for a charmless B meson 
decay that addresses both the construction of the operator basis and the factorization proper¬ 
ties of the matrix elements in a systematic way. The same techniques can be used to establish 
factorization for hadronic decays such as B ^ tttt. 
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A Matching coefficients 


In the matching onto the SGETj current operators, the following combinations of Wilson 
coefficients appear: 


^ /~iA c^A _ 


_41,i"-in_ioinA_4ln^:SSE 

xm^ rrib xrrib 


2 + 


ttsCf 

471 


4Li2(l -x)-12 






2x + 0{as ), 
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leg + eg = 0{a .), 


where x = 2E/mb- The A-type coefficients can be found in [2, 33]. The 0{as) term for the 
S-type coefficients are also known, see [33, 34], However, they appear only at 0{al) in the 
decay rate, since they are multiplied by the jet function, which is proportional to 

We collect the known results from the literature for the matching coefficients appearing 
in Section 4. For the matching of Qi and Qs onto A-type current operators, we need the 
functions 


Gi{x) 


8 1 A'-qcd 11 27ri 27r^ 

- In H-^- 

3 mb 3 3 9 

104 WocD 833 20f7r o/o 

-In -^- 

27 mb 162 27 9 

2 r 

- 48 -|- 30f7r — — 2m^ — 36^(3) -|- (36 -|- — 97r^) In a; 


-I- (3 -I- 6f7r) \v?x -I- In^x 


X 


+ 

+ 


2 

9 

1 


18 -|- 27r^ — 2m^ -|- (l2 — 67r^) Inx -|- 6f7rln^a; -|- In^x 
— 9 -|- 112i7r — 147r^ -|- (182 — 48f7r) In a; — 126 In^a; 


X 


2 


X^ + 0{x‘^) , 
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which were dehned in [18] and deduced from the results of [35, 36]. Here a; is a ratio of 
squared masses, Xq = my ml, for a quark of flavor q. The matching of Qi onto B- and G-type 
operators involves the functions 


1 -I- 4a; 


fix) = 


.vr 


l2 


arctanh(Vl — 4a;) — i— \ ; for a; <1/4, 


1 — 4a; arctan 


V4a; — 1 ’ 


for X > 1 /4 , 


and [38] 


G(a;, m) = —4 / dvv{l — v)\n[x — v{l — v)u\ 
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